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Abstract.  We  examine  one  way  in  which  a  compliant  boundary  may  affect  the 
low  wave  number  pressure  field  measured  at  the  surface  under  a  turbulent 
boundary  layer.  Specifically,  we  assume  that  the  low  wavenumber  pressure 
field  measured  at  the  wall  results  from  the  large  scale  motions  in  the 
outer  part  of  the  flow,  and  that  these  are  triggered  by  the  bursting  in  the 
buffer  layer,  which  is  caused  in  its  turn  by  secondary  instabilities  grow¬ 
ing  on  inflectionary  profiles  caused  by  longitudinal  vortices  in  the  sub- 
and  buffer  layers.  We  examine  the  direct  effect  of  the  compliant  boundary 
on  the  growth  of  these  longitudinal  vortices,  using  the  energy  method.  We 
represent  the  compliant  boundary  by  a  half  space  filled  with  a  linear 
visco-elastic  medium  with  a  single  time  constant.  We  find  that  the  effect 
of  the  boundary  is  to  introduce  an  effective  slip  in  the  boundary  condition 
on  the  disturbance  (the  longitudinal  vortices),  producing  a  boundary  condi¬ 
tions  intermediate  between  the  inviscid  boundary  condition  and  the  no-slip 
boundary  condition.  Since  this  boundary  condition  is  less  dissipative,  the 
vortices  are  less  stable,  and  have  a  higher  growth  rate.  They  might 
consequently  be  expected  to  grow  faster  and  occur  more  frequently,  if  this 
were  the  only  mechanism  involved,  resulting  in  a  shift  of  the  low  wave- 
number  pressure  spectrum  to  higher  frequencies.  This  appears  to  be  the 
opposite  of  observations,  suggesting  that  other  mechanisms  are  involved. 


1.  Scientific  background.  As  early  as  Lumley  (1964)  the  possibility  was 
considered  that  a  compliant  surface  might  alter  the  structure  of  a  turbu¬ 
lent  boundary  layer.  The  interest  at  that  time  was  in  drag  reduction;  it 
was  hoped  that  a  compliant  surface  might  change  the  stability  characteris¬ 
tics  of  the  large  eddies  in  the  viscous  sublayer  and  buffer  region,  result¬ 
ing  in  a  thicker  sublayer,  corresponding  to  a  reduction  of  drag.  A  prelim¬ 
inary  analysis  was  carried  out,  which  indicated  that  realistic  compliant 
surfaces  were  probably  always  detrimental.  The  analysis,  however,  was 
limited  to  progressive  waves;  a  treatment  of  the  more  realistic  longitudi¬ 
nal  vortices  had  been  possible  since  Serrin  (1959),  but  Lumley  was  unaware 
of  this  work  until  1967  (see  Lumley,  1971).  At  the  time  (1964)  the  matter 
was  not  pursued,  since  other  forms  of  drag  reduction  appeared  more  inter¬ 
esting.  The  subject  has  recently  been  reexamined  (Bushnell  et  al,  1977), 
and  considerable  experimental  work  exists  (some  of  it  indicating  reduction 
of  drag);  the  theoretical  picture  is  not  clear,  although  a  modification  of 
the  large  eddy  flow  preceding  a  burst  is  thought  to  be  implicated. 

A  possibility  that  was  overlooked  at  the  time,  but  that  may  be  of 
interest  now,  is  that  alteration  of  the  structure  of  the  turbulent  boundary 
layer  by  a  compliant  surface  might  alter  the  pressure  field  observed  at  the 
surface,  particularly  in  the  low  wavenumber  region,  i.e.,  the  region  of 
disturbances  having  scales  of  many  boundary  layer  thicknesses.  This  would 
permit  the  modification,  by  the  application  of  a  compliant  treatment  to  a 
rigid  surface,  of  the  low  wavenumber  pressure  field  perceived  at  the  rigid 
surface,  something  not  otherwise  possible  without  structural  modification, 
since  the  application  of  a  relatively  thin  compliant  surface  treatment 
would  not  influence  the  transmission  of  pressure  fields  in  this  wavenumber 
range. 

Examination  of  such  a  possibility  can  more  easily  be  undertaken  now 
for  several  reasons:  first,  the  results  of  Serrin  (1959),  and  in  particular 
Lumley  (1971),  already  referred  to,  which  make  possible  the  treatment  of 
the  longitudinal  vortices;  second,  an  understanding  of  the  dynamical  role 
of  these  vortices,  arising  in  part  from  studies  of  polymer  drag  reduction 
(Lumley,  1977);  and  finally,  the  existence  of  much  recent  data  on  wall 
pressure  fluctuations  and  their  relation  to  boundary  layer  structure  (see 
Willmarth  &  Bogar,  1977). 
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The  mechanism  we  have  in  mind,  to  connect  the  large  eddies  in  the  sub- 
and  buffer  layers  with  the  pressure  distribution  on  the  wall  is  this  (see 
Willmarth  &  Bogar,  1977):  the  large  eddies  are  in  the  form  of  counter¬ 
rotating  streamwise  vortices;  in  the  updraft  between  two  vortices,  the 
streamwise  velocity  profile  develops  an  inflection  point,  due  to  the  slow 
moving  fluid  from  the  vicinity  of  the  wall  which  is  lifted  up  into  the 
higher  velocity  region;  the  inviscid  instability  resulting  from  this 
inflectionary  profile  produces  a  burst  of  Reynolds  stress;  this  burst  of 
Reynolds  stress  results  in  a  temporary  change  in  the  thickness  of  the 
boundary  layer,  associated  with  a  fluctuation  in  the  pressure  which  is  felt 
throughout  the  thickness  of  the  layer,  and  in  particular  at  the  wall. 

Hence,  the  key  to  the  whole  mechanism  is  the  frequency  of  bursting  and 
the  growth  rate  of  the  large  eddies  in  the  sub-  and  buffer  layers.  Exist¬ 
ing  data  is  reasonably  consistent  regarding  the  scaling  of  the  size  of 
these  large  eddies,  which  are  observed  to  scale  with  the  wall  variables. 
Unfortunately,  the  scaling  of  the  frequency  of  bursting  is  less  secure: 
there  is  data  to  indicate  scaling  both  with  wall  variables  and  with  outer 
layer  variables.  It  is  thus  not  at  all  clear  what  mechanism  controls  the 
bursting  frequency.  We  will  examine  here  the  possibility  that  the  growth 
rate  of  the  eddies  controls  the  bursting  frequency;  that  is,  that  the 
life-time  of  the  eddy  is  inversely  proportional  to  the  growth  rate  during 
the  time  of  most  rapid  growth;  that  the  eddies  are  amplified  from  a  primor¬ 
dial  soup  of  disturbances  in  the  sublayer  (containing  all  possibilities) 
and  that  consequently  a  new  eddy  will  appear  as  soon  as  possible  after  the 
last  one  has  died;  hence,  that  the  frequency  of  appearance  will  be 
controlled  by  the  life-time. 

In  order  to  determine  the  growth  rate  of  the  eddies,  we  will  use  the 
energy  method  (Lumley,  1971).  In  this  method,  we  ask  for  the  form  of  the 
disturbance  which  has  the  largest  growth  rate  at  a  given  instant,  subject 
to  given  boundary  conditions,  and  for  a  given  mean  velocity  profile.  Such 
an  optimal  disturbance  will  not  remain  optimal  as  it  evolves,  and  was  not 
optimal  before  the  instant  of  analysis.  Among  other  effects,  the  presence 
of  the  disturbance  will  progressively  erode  the  mean  velocity  profile  (due 
to  the  Reynolds  stress  produced  by  the  disturbance),  and  the  reduced  strain 
rate  of  the  mean  profile  will  reduce  the  energy  supply  of  the  disturbance, 
slowing  its  growth  rate,  and  increasing  its  size.  We  will  not  analyze  this 
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non-linear  effect  here,  which  has  not  yet  been  analysed  for  the  rigid  wall, 
but  will  content  ourselves  with  an  analysis  of  the  influence  of  the  bound¬ 
ary  condition  corresponding  to  a  compliant  surface  on  the  growth  rate  of 
these  optimum  disturbances,  so  that  we  can  have  a  direct  comparison  with 
Lumley  (1971). 

2.  Extremization  with  variable  viscosity.  We  represent  the  increased 
momentum  transport  of  the  turbulence  in  the  boundary  layer  by  a  variable 
viscosity;  this  has  known  drawbacks  (Tennekes  &  Lumley,  1972),  but  Lumley 
(1971)  shows  that  very  little  of  the  turbulent  momentum  transport  enters 
into  the  dynamics.  Specifically,  the  dynamics  are  dominated  by  the  region 
very  near  the  wall,  where  the  strain  rate  and  the  momentum  transport  are 
large,  and  the  latter  is  predominantly  viscous. 

We  begin  from  the  Navier-Stokes  equations  written  with  variable 
viscosity,  for  the  undisturbed  motion  U£  and  for  the  combination  of 
undisturbed  motion  plus  disturbance,  u£.  Multiplying  by  the  disturbance 
velocity  and  summing,  an  equation  for  the  disturbance  energy  e  ■■  u^u^/2 
can  be  written.  This  equation  is  now  integrated  over  the  region  in  which 
the  flow  is  taking  place;  we  assume  that  the  undisturbed  boundaries  are  not 
moving  (at  least  not  normal  to  themselves).  Finally,  we  introduce  the 
notation 


o/edV  -  (d/dt)/edV 


(2.1) 


Using  the  divergence  theorem  to  rewrite  all  transport  terms  as  surface 
integrals,  we  finally  have 


o/edV 


■  “/s. .u.u.dV  -  2/vs. .s 

^  ij  1  j  ij 

-  ^[u.e  -  u.T. ./pjdo. 
J  1  J 


.  .dV 


(2.2) 
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where  S.  .  is  the  mean  strain  rate  (U.  .  U .  .  )/2  and  s..  the  strain 
ij  i.J  J.i 

rate  of  the  disturbance  motion  (u.  .  +  u.  .)/2.  t..  ■  -p6.  .  +  2ps.. 

i.J  J.i  ^ 

is  the  stress  tensor  of  the  disturbance  motion. 

The  volume  integral  may  be  over  the  entire  region  in  which  the  flow  is 
taking  place.  However,  if  the  region  is  unbounded  in  one  or  more  direc¬ 
tions,  we  will  suppose  that  the  motidn  is  periodic  in  that  direction,  and 

we  will  take  the  integral  only  over  a  period. 

The  first  term  in  the  surface  integral  represents  the  transport  of 
disturbance  energy  across  the  boundary  by  the  disturbance  velocity.  The 
contribution  to  the  integral  from  the  parts  of  the  surface  that  lie  in  the 
flow,  in  the  case  of  a  periodic  motion,  and  that  are  consequently  at  the 
beginning  and  end  of  a  period,  will  cancel,  so  that  we  need  concern  our¬ 
selves  only  with  that  part  of  the  integral  which  lies  on  an  actual  bounding 
surface.  We  can  eliminate  this  term  using  either  of  two  arguments.  First, 
we  suppose  that  the  surface  is  nearly  rigid,  so  that  the  motion  of  the 
surface  is  relatively  small.  Essentially,  we  can  imagine  expanding  in  a 
parameter  related  to  the  surface  rigidity.  Then  the  disturbance  velocity 
at  the  surface  is  a  small  quantity,  which  vanishes  when  the  surface  is 
rigid.  The  surface  stress  is  of  order  unity,  since  it  does  not  vanish  as 
the  surface  becomes  rigid.  Consequently,  the  second  term  in  the  surface 

integral  is  of  first  order,  while  the  first  terra  is  of  third  order,  and  can 

be  neglected  by  comparison. 

A  second  argument,  which  is  unrelated  to  the  amplitude  of  the  surface 
motions,  supposes  (as  we  will  suppose)  that  the  disturbance  is  periodic 
cross-stream.  Then  the  disturbance  energy,  being  quadratic,  will  contain 
components  at  zero  wavenumber  and  at  twice  the  wavenumber  of  the  basic 
disturbance;  consequently,  the  energy  flux  which  appears  in  the  surface 
integral  will  contain  components  at  the  basic  wavenumber  and  at  three  times 
the  wavenumber,  all  of  which  will  integrate  to  zero  over  one  or  more 
periods  of  the  basic  disturbance. 

We  must  now  extremize  equation  (2.2)  (suppressing  the  energy  flux  term 
in  the  surface  integral).  We  wish  to  consider  only  incompressible  motions, 
so  we  must  add  a  Lagrange  multiplier  to  the  equation  to  maintain  incom¬ 
pressibility,  which  will  take  the  form  of  a  pressure.  We  add  to  the  right 


hand  side  the  following; 
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0  -  2/{(()/p)'5u.  .dV  -  2^(<ti/p)6u.da.  -  2/((j)  ./p)6u.dV  (2.3) 


where  t  is  the  Lagrange  multiplier.  This  gives  for  the  extremized  form 


/(lu.  +  2S..U.  -  4(vs..)  .  +  2(j>  ./p)6u.dV 
1  J  .J  .1  1- 


^(6u.t.,/p  u.6t../p  -  2Tt,fiu./p)do. 

1.  ij  i  tj  ij  1  j 


(2.4) 


where 


xt,  »  ,  +  2U8..  (2.5) 

ij  ij 


We  may  presume  that  the  variation  in  the  disturbance  velocity  in  the 
interior  of  the  fluid,  and  at  the  surface,  are  independent  (since  the 
variations  at  any  two  points  are  independent);  hence,  the  volume  integral 
gives  (under  the  usual  assumption  that  the  integrand  is  continuous) 


ou. /2  +  S. .u. 
J-  J 


-<t>  j/p  + 
$  X 


2( vs. .) 


(2.6) 


We  cannot  remove  the  integrand  from  the  surface  integral  yet,  because  we 
have  not  related  the  variation  in  the  surface  stress  to  the  variation  in 
the  disturbance  velocity. 
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3.  The  Viscoelastic  Substrate.  We  begin  by  assuming  that  the  wall  is  a 
semi-infinite  region,  a  half-space.  There  are,  of  course,  other  possibili¬ 
ties:  we  could  consider  a  visco-elastic  slab  of  finite  thickness  over  a 
rigid  surface,  or  other  wall  compositions  with  complex  compliance  varying 
with  depth  in  different  ways.  The  principal  way  in  which  such  walls  would 
differ  from  the  half-space  we  have  chosen,  would  be  in  the  existence  of 
reflection  from  the  compliance  changes.  Non-dimensionalized  with  sublayer 
variables,  the  time  required  for  a  disturbance  to  travel  through  the  sub¬ 
strate  a  distance  equal  to  the  sublayer  thickness  is  about  11.6,  if  Che 
wall  material  is  of  the  same  density  as  Che  fluid,  and  Young's  modulus  is 
equal  Co  Che  dynamic  head  based  on  Che  shear  velocity.  Non-dimensionalized 
in  Che  same  way,  the  time  required  for  substantial  growth  of  Che  eddy  is 
roughly  4,  using  values  corresponding  to  a  rigid  wall.  This  is,  of  course, 
a  very  compliant  wall;  as  the  wall  is  stiffened,  Che  transmission  time  will 
be  reduced  (proportional  to  the  inverse  of  the  square  root  of  Young's 
modulus).  At  Che  same  time,  as  Che  eddy  grows,  and  Che  shear  is  reduced, 
the  growth  rate  will  decrease,  increasing  the  growth  time.  Hence,  the 
times  are  probably  comparable,  with  the  reflection  time  likely  being  Che 
shorter  for  stiffer  walls.  Thus,  there  is  a  possibility  that  reflection 
plays  a  role,  and  this  bears  looking  into.  Here,  however,  we  concentrate 
on  Che  mechanisms  present  without  reflection.  We  will  find  chat  Che  prin¬ 
cipal  mechanism  has  to  do  with  Che  displacement  of  Che  surface,  which 
requires  a  compensating  disturbance  velocity  Co  cancel  the  undisturbed 
velocity  at  Che  new  location  of  Che  surface  (see  section  4.) 

If  we  write  Sj  for  the  diplacement  vector  at  each  point  in  the 
substrate,  Che  equations  of  motion  may  be  written  as 


S.  =  -p  ./p’  +  T. .  ./p’ 


(3.1) 


»  V  *• 

Primed  parameters  refer  to  wall  materials.  The  stress  tensor  in  the 

substrate  is  given  by 
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tY.  -  G(S.  .  +  S.  .)/2  +  w'(d/dt)(S.  .  +  S.  .)/2  (3.2) 

ij  i.J  J.i  i.J  J.i 


In  our  extremization  procedure  we  say  nothing  about  the  temporal 
evolution  of  the  disturbance.  The  quantity  a  which  appears  in  the  equa¬ 
tions  is  the  instantaneous  relative  growth  rate  of  the  integrated  energy. 
There  is  no  implication  that  the  growth  rate  is  the  same  at  different 
places  in  the  fluid,  or  will  be  the  same  at  other  instants  before  or  after 
the  instant  for  which  the  equations  are  written.  In  fact,  Lumley  (1971) 
shows  that  Che  solution  to  Che  energy  method  equations  necessarily  has  a 
non-linear  evolution,  with  different  growth  rates  at  different  points,  and 
that  it  is  optimum  only  at  the  instant  of  analysis,  evolving  through 
non-optimum  states  before  and  after,  and  hence  necessarily  having  different 
growth  rates.  In  order  Co  obtain  the  stresses  in  the  substrate,  however, 
we  must  have  a  value  for  the  relative  growth  rate  of  the  displacement,  and 
f-;r  its  second  derivative.  In  the  absence  of  other  information,  the 
simplest  choice  appears  to  be  to  suppose  that  the  displacements  in  the 
substrate  are  growing  exponentially  in  a  neighborhood  of  the  instant  of 
analysis,  and  that  the  relative  growth  rate  is  one  half  of  the  relative 
growth  rate  of  the  integrated  disturbance  energy  (the  same  as  the  growth 
rate  of  the  disturbance  velocity,  if  it  were  everywhere  the  same).  Hence, 
we  write 


S.  =  S^a/2,  »  Sj^(a/2)^ 


(3.3) 


which  perm..ts  us  to  write 


tY.  ”  (G  +  u'a/2)(S.  .  +  S.  .)/2 
J-J  i.J  J.i- 


(3.4) 
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As  a  disturbance  in  the  fluid,  we  will  be  interested  in  a  velocity 
field  that  is  not  a  function  of  the  streamwise  coordinate,  and  is  periodic 
cross-stream,  because  this  is  the  type  of  disturbance  which  is  observed  in 
the  sublayer  over  a  rigid  wall  (see  Lumley,  1971)  as  well  as  over  a 
compliant  one  (see  Hansen  et  al,  1980).  Consequently,  we  will  be  inter¬ 
ested  in  disturbances  of  the  substrate  of  the  same  type,  and  we  will 
assume 


S 


2 


0 


f  *  ^exptikx^] , 


A 

Sj  =  Sexp[ikX2] 


(3.5) 


Taking  account  of  the  radiation  condition  (requiring  that  the  disturb¬ 
ances  in  the  substrate  vanish  at  infinity),  the  solutions  to  the  equations 
of  motion  may  be  written  as  (in  terms  of  arbitrary  boundary  conditions) 

?  »  (X$Q  -  yk)exp[kX2]/(X  -  1)  +  (yk  -  ¥Q)exp[kXx2]/(X  -  1) 

(3.6) 

S  *  SQexp[kXx2l 


where 


X  -  [1  +  a2p'/2k^(G  +  u’a/2)]^^^ 


(3.7) 


The  stress  tensor  at  the  surface  may  now  be  written  in  terms  of  the  arbi¬ 
trary  boundary  conditions,  where  we  obtain  the  pressure  from  the  cross 
stream  component  of  (3.1): 


■MMKi 


-Pq  -  i(G  +  w'o/2)(X'Fq  -  +  X)/2 


(3.8) 


so  that  the  symmetric  stress  is  given  by  t.  ,  =  -S^'S.  .  +  * 

ij  ij  ij 

id  +  X)(X$Q  -  ?^/k)k2/2  •  •  • 

XSk/2  ‘  iX(l  +  X)?Qk2/2  +  id  -  X)$^k/2  *  *  * 

iSk/2  {(X  -  Dfg  -  (1  +  X)^yk}k2/2  iXCl  +  X)fQk2/2  -  i(3  +  X)?^k/2 


X  (G  +  u'o/2) 


(3.9) 


Now,  we  are  interested  in  the  surface  response,  so  that  we  may  write 
on  the  surface 


T.  .do.  *  g.,  S.do 
ij  J  Ik  k 


(3.10) 


which  becomes,  on  the  bottom  surface 


^2  ’  *ik\®^^ 


(3.11) 


From  (3.9)  and  the  definition  (3.7),  then,  we  may  write 


X/(X2  -  1) 


V(X-l)  i/(X+l)  )  p'a/2k  (3.12) 


-i/(X  +  1)  1/(X  -  1) 
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4.  The  boundary  condition.  Since  the  surface  of  the  viscoelastic  substrate 
moves,  it  interacts  with  the  fluid  motion  other  than  simply  through  the 
stress  on  the  surface.  In  particular,  suppose  that  the  surface  moves 
normal  to  itself  into  the  fluid,  without  moving  laterally.  The  no  slip 
condition  must  still  be  maintained;  the  undisturbed  velocity  at  the  new 
location  of  the  surface  in  now  non-zero  and  positive,  so  that  a  negative 
streamwise  disturbance  velocity  at  this  location  is  required  to  cancel  the 
former. 

The  classical  kinematic  boundary  condition  takes  a  somewhat  different 
form  here;  the  part  of  the  condition  relating  to  the  slope  of  the  surface 
is  second  order,  and  may  be  neglected;  at  the  same  time,  we  have  a  new 
condition  resulting  from  the  possibility  of  lateral  motion.  The  general 
condition  becomes: 


§i  -  U.(Xj  +  Sj,  82,83)  =  6j,U(S2)  +  u.(xj  +  Sj,  $2,  S3) 


(4.1) 


.U'S-  +  u?  +  u?  .8.  + 
li  2  1  i,j  j 


Now,  if  the  wall  were  rigid,  the  disturbance  velocities  at  the  surface 
would  vanish,  as  would  their  gradients  in  the  plane  of  the  surface;  hence, 
most  of  the  third  and  subsequent  terms  in  (4.1)  will  be  of  second  or  higher 
order,  and  may  be  neglected  if  the  surface  displacements  are  small.  The 
gradients  normal  to  the  wall  of  the  disturbance  velocities  parallel  to  the 
wall  do  not  vanish  as  the  wall  displacement  vanishes,  however,  and  hence 
may  not  be  taken  to  be  small.  To  avoid  the  introduction  of  non-linearities 
(which  do  not  seem  to  be  worth  the  additional  effort  they  would  require 
since  we  expect  these  results  to  be  of  at  most  qualitative  significance), 
we  will  assume  that  the  disturbance  shear  is  small  relative  to  the  mean 
shear.  Then  the  relation  between  the  surface  displacements  and  the 
disturbance  velocities  at  the  surface  becomes: 

Sjo/2  -  i322U’/a  +  ilj 
S2O/2  -  d® 

S3  0/2  -  u° 


(4.2) 


Note  that  the  lateral  surface  displacement  may  not  be  in  the  expected 
direction,  corresponding  to  the  value  of  the  lateral  disturbance  velocity 
in  the  fluid  at  the  surface,  if  the  vertical  disturbance  velocity  (combined 
with  the  undisturbed  shear  and  the  growth  rate)  has  the  right  value.  We 
will  find,  in  fact,  that  this  is  a  dominant  effect. 

Now,  on  the  bottom  surface  we  may  write 


g.j(J°2UVo  +  +  g.2«i2  +  gi3u^3  » 


(4.3) 


We  may  now  return  to  the  surface  integral  part  of  (2.4).  If  we 
suppose  (corresponding  to  the  type  of  disturbance  in  which  we  are  inter¬ 
ested)  that 


Use  of  (3.12)  and  (4.3)  permits  evaluation  of  the  coefficient  matrix  to 
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X/(X2  -  1) 
XU'/o(x2  - 

0 


1) 


XU'/o(x2  -  1) 

X/(X  -  1) 

-i/(X  +  1) 


0 

i/(  X  +  1)  \  p'  a/2k 
1/(X  -  1)  / 

(4.7) 


5.  Solution  of  the  stability  problem.  We  suppose  a  disturbance  of  the  form 

Uj  *  ilexp[ikx^] 

U2  =  'll  2  *  ik({iexptikx2l  (5.1) 

U2  ”  -'ll  2  *  -IJi'expfikx^] 

We  may  obtain  an  expression  for  the  pressure  at  the  surface  from  the 
lateral  component  of  equation  (2.6),  giving 


-f  «  (-t|i'a/2  +  V  (D^  -  k^)(|i']p/ik  (5.2) 

m 

where  a  subscript  zero  indicates  evaluations  at  the  wall. 

Using  this,  (5.1)  and  (4.7),  we  may  evaluate  the  stress  in  the  fluid,  and 
obtain  the  boundary  conditions  explicitly  as 

«  [i3°X/(X^  -  1)  +  ikfp^U' X/a(X^  -  l)]p'o/2k 

l-f'a/2  +  \^(D^  -  3k^)ip°']p/ik  =  ['3°XU'/o(X^  -  1)  +  ikip°X/(X  -  1) 

-i$0'/(X  l)]p'o/2k 

-p(D^  ♦  k^)(fi®  =  [k<f°/(X  +  1)  -  f'/iX  -  l)]p’a/2k 


We  wish  now  to  introduce  the  same  non-dimensionalization  as  used  in 
Lumley  (1971);  that  is,  we  non-dimens ionalize  lengths  by  the  thickness  of 
the  viscous  sublayer  (defined  as  the  point  of  intersection  between  the 
sublayer  profile  and  the  logarithmic  profile),  and  introduce  as  a  small 
parameter  c  the  inverse  of  the  Reynolds  number  defined  using  this  length, 
the  shear  velocity  and  the  wall  viscosity.  The  dimensionless  growth  rate 
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is  defined  as  3^  >  o/U' .  Indicating  dimensionless  variables  by  the  same 
letters,  the  boundary  conditions  become 

2ke^a°’  -  3^(a°  +  -  1) 

e^(D^  -  3k^)(|>°’  =  0^iiXG°/8^(X^  -  1)  -  k?»°X/(X  -  1) 

+  ip^'d  +  (X  +  l)p/p']/(X  +  l)}p'/2p  (5.4) 

-2e^(D^  +  k^)(I.°  -  6^{k?»°/(X  +  1)  -  (P°V(X  -  l)}p'/p2k 
The  expression  for  X  may  be  written  as 


X^  =  1  +  8^p7p2k^e^(G/pu^  +  g^p'v'/pv  ) 

"  ni 


(5.5) 


where  primed  parameters  correspond  to  wall  material.  We  have  in  mind 
constructing  an  asymptotic  expansion  for  large  Reynolds  number,  so  that  the 
region  of  viscous  influence  on  the  disturbance  motion  will  become  thin  and 
can  be  handled  by  matching,  as  in  Lumley  (1971).  There  are  a  number  of 
ways  that  the  compliant  wall  can  be  incorporated  in  this  scheme,  corre¬ 
sponding  to  different  limiting  processes.  Since  in  any  event  we  are  inter¬ 
ested  primarily  in  learning  the  qualitative  effects  of  the  compliance,  we 
would  be  satisfied  by  some  sort  of  expansion  about  the  rigid  wall  case. 
The  easiest  way  to  accomplish  this  is  to  take  the  Young's  modulus  propor¬ 
tional  to  the  square  of  the  Reynolds  number,  so  that  the  wall  becomes  more 
and  more  rigid  as  the  Reynolds  number  increases.  This  has  the  effect  of 
holding  X  constant  as  the  Reynolds  number  goes  to  infinity.  Then  if  we 
write 


u  ■  ^^(hQ  ♦  e^hj  ♦  . . . ) 


(5.6) 


we  find  that  the  lowest  order  solution  corresponds  to  the  rigid  wall  case, 
with  ■  gg  ■  gg'  ■  0  at  the  wall.  To  second  order  we  obtain 
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8^(h°  +  g°ik/3^)Xp'/p(X^  -  1)  -  2kh®'  - 

e^[ih°X/8^(X^  -  1)  -  kg®X/(X  -  1)  +  g°’(l  +  (X  +  l)p/p’] 

A  .V  ^  -  (5.7) 

/(X  +  l)]pV2p 

B^[kg°/(X  +  1)  -  gJ'/(X  -  l)]p'/p2k  -  -2gJ"  -  2DB^ 

The  expressions  for  the  lowest  order  solutions  will  now  be  the  same  as  in 
the  rigid  wall  case,  and  may  be  taken  directly  from  Lumley  (1971).  These 
have  been  used  to  give  the  final  expressions  in  (5.7),  where  D  is  a  nega¬ 
tive  constant  which  must  be  determined  by  matching. 

We  may  solve  (5.7)  to  obtain 

h°  -  -2iD(X  -  1){1  +  k^(X  +  l)/40^  -  2k(X  -  1)[1  +  (X  +  l)p/p']/8(X  ♦  l)}p 
/p'6{X[l/3^(X  ♦  1)  -  1]  +  (X  -  1)2[1  +  (X  ♦  l)p/p’]/(X  +  1)2}  (5  gj 

g°  *  D(X  -  l){2e  +  k^/2B^  +  4k(X  -  1)[1  +  (X  +  l)p/p']/(X  +  l)}p 

/p'k{X/0^(X  +1)  -  X  +  (X  -  1)2[1  +  (X  +  l)p/p’]/(X  +  1)2}  (5.9) 

gj'  =•  -D(X  -  l){4k  -  (X  -  1){2B  +  k^/23^  +  4k(X  -  1) 

[1  +  (X  +  l)p/p']/(X  +  1)}  /(X  +  1){X/3^(X  +  1)  -  X  (5.10) 

+  (X  -  1)^[1  +  (X  +  l)p/p’]/(X  +  l)^}}p/p' 

In  Lumley  (1971)  it  is  shown  that,  in  order  to  have  one  turning  point 
in  the  outer  solution,  so  that  the  eddies  close  toward  the  outside  of  the 
layer  (rather  than  producing  upward  and  downward  jets  which  persist  to  the 
outside  of  the  layer),  it  is  necessary  to  have  1/23^  >  2.  With  this 
restriction,  it  is  clear  from  (5.8-10)  that  g°'  is  positive  real,  g°  is 
negative  real  and  h°  is  positive  imaginary.  (We  presume  that  the  outer 
stream  function  is  positive  real,  corresponding  to  D  negative).  In 
figure  1  we  indicate  the  inviscid  (outer)  stream  function  and  its  first 
derivatives,  together  with  the  amplitude  of  the  streamwise  disturbance;  we 
have  shown  dotted  the  modification  due  to  viscosity  over  a  rigid  wall,  and 
the  dashed  curve  indicates  the  modification  due  to  compliance.  It  is  seen 
that  the  compliant  solution  is  intermediate  between  the  inviscid  and  the 
viscous  solution.  We  would  thus  expect  the  compliant  solution  to  be  less 


dissipative  (more  like  the  inviscid  solution)  and  hence  more  unstable. 
Note  that,  although  the  lateral  velocity  at  the  wall  is  in  the  direction 
that  might  be  expected  based  on  the  outer  solution,  both  the  normal  veloc¬ 
ity  and  the  streamwise  velocity  are  reversed.  This  may  be  attributed 
directly  to  the  effect  of  the  boundary  condition  in  section  4. 

Now,  we  must  obtain  the  solution  for  the  g^  corresponding  to  the 
boundary  conditions  (5.8-10)  (since  we  will  do  the  matching  to  obtain  the 
stability  criterion  entirely  in  terms  of  gj ,  we  will  not  need  hj). 
Since  the  rigid-wall  boundary  conditions  on  g]^  vanish,  it  is  convenient 
to  write  g^  as  the  sum  of  the  rigid  wall  solution  and  a  particular  solu¬ 
tion  corresponding  to  the  non-zero  boundary  conditions,  which  is  easily 
shown  to  be 


*1 


(5.11) 


(where  the  ellipsis  indicates  the  rigid-wall  solution).  The  matching  may 
now  be  carried  out  exactly  as  in  Lumley  (1971),  resulting  in  an  additional 
term  on  the  left  of  the  O(e^); 


-B  +  g 


(5.12) 


(the  ellipsis  indicating  the  unaltered  right  hand  side  of  the  corresponding 
equation  in  Lumley,  1971)  and  on  the  left  of  the  0(e^5): 


8[B  +  (9/4  +  4B^)k^D/246^1  +  g°’  -  ... 


(5.13) 


Adding  these  expressions  results  in  one  additional  term  on  the  right  hand 
side  of  U3: 


“3 


. ..  -  k(l/4B^  -  l)^^^(gj  +  gJV8)/DB 


(5.14) 
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Finally,  the  eigenvalue  condition  is  obtained; 


3ir/4  -  Y_Qk 


A  1/2  A  1/2  A  .33,3 

-(1/46^  -  1)  kc/6  -  (1/46^  -  1)  (1/486^  +  Ok'^e-’/ze 

-  (1/48^  -  l)^^^gJ  +  gJ'/6)e\/De 


This  condition  is  the  same  as  in  Lumley  (1971)  except  for  the  additional 
term  on  the  right  hand  side. 

For  a  fixed  value  of  the  growth  rate,  the  condition  (5.15)  gives  in 
general  two  values  of  wavenumber  for  each  value  of  Reynolds  number  (1/e); 
there  is  a  minimum  value  of  Reynolds  number,  corresponding  to  a  critical 
value  of  wavenumber.  A  typical  curve  is  shown  in  figure  2.  In  Lumley 
(1971)  this  minimum  value  was  found  by  differentiating  (5.15)  with  respect 
to  wavenumber,  and  requiring  that  both  the  value  and  the  derivative  vanish. 
Here,  however,  the  situation  is  algebraically  so  complex  (due  to  the 
appearance  of  wavenumber  in  the  expressions  for  surface  response)  that  it 
was  felt  to  be  simpler  to  determine  the  minimum  value  of  Reynolds  number 
numerically  from  (5.15)  for  each  fixed  value  of  growth  rate.  The  resulting 
plot  is  shown  in  figure  3,  where  the  rigid  wall  case  is  also  shown.  It  is 
seen  that  the  growth  rate  curve  is  shifted  upward,  so  that  a  larger  growth 
rate  is  experienced  for  a  given  Reynolds  number  with  a  compliant  surface 
than  with  a  rigid  surface,  as  was  anticipated  from  the  modification  of  the 
boundary  condition,  which  was  seen  to  be  less  dissipative. 

6.  Discussion.  Lumley  (1971)  is  included  as  an  appendix  to  this  report  for 
convenience.  Comparison  of  the  rigid  wall  curve  in  figure  3  with  the  cor¬ 
responding  curve  in  Lumley  (1971)  indicates  a  considerable  difference, 
which  has  been  traced  to  a  numerical  error  in  the  computations  made  for  the 
curve  in  Lumley  (1971).  The  equations  in  Lumley  (1971)  produce  the  curve 
in  figure  3,  which  is  correct.  The  curve  in  figure  3  may  be  rationalized 
physically  in  the  following  way:  note  that  the  wavenumber  is  roughly  the 
same  as  that  for  the  Couette  flow  (the  left-hand  curve  in  Lumley  1971, 
which  is  correct);  hence,  under  the  same  shear,  the  eddies  are  about  the 
same  size,  whether  an  upper  wall  is  present  or  not.  The  Reynolds  number, 
on  the  other  hand,  is  considerably  increased,  though  not  doubled;  this  is 
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because  Che  outer  half  of  Che  eddy  is  unsupported  energetically  -  the  shear 
in  this  region  is  very  much  reduced,  so  Chat  this  part  of  the  eddy  must  be 
supported  dynamically  by  Che  part  near  Che  wall. 

Note  Chat  Che  tendency  of  Che  curves  in  figure  3  to  bend  back  Co  the 
right  as  they  approach  the  abscissa  is  not  real,  but  results  from  the 
failure  of  Che  various  approximations  (high  wave  number,  high  growth  rate) 
Chat  were  made  in  Che  course  of  the  calculation.  The  curves  should  be 
extrapolated  in  Che  same  way  as  were  the  curves  in  Lumley  (1971). 

The  lambda  indicated  on  figure  3  is  the  same  as  that  used  in  Lumley 
(1971),  and  is  not  the  lambda  of  this  report.  The  lambda  of  figure  3  is 
Che  growth  rate  non-dimens ionalized  by  sublayer  variables. 

The  very  small  difference  of  figure  3  between  the  rigid  and  compliant 
wall  (despite  Che  very  compliant  values  of  Che  parameters  chosen)  is  disap¬ 
pointing,  as  is  Che  fact  that  it  is  in  Che  wrong  direction.  We  must  con¬ 
clude  Chat  substantial  interaction  with  a  compliant  surface  does  not  exist 
with  eddies  of  this  form.  Probably  only  for  a  progressive  wave  motion  does 
a  strong  interaction  exist.  The  fact  that  the  change  is  in  the  wrong  di¬ 
rection  results  from  the  fact  Chat  the  primary  effect  of  the  compliant 
surface  here  is  to  allow  an  apparent  slip  at  Che  wall,  due  to  the  interac¬ 
tion  of  the  surface  motion  and  the  undisturbed  shear.  If  we  are  to  have  a 
favorable  change,  we  must  seek  another  mechanism. 

Several  possibilities  come  Co  mind.  First  is  the  possibility  that  Che 
compliant  surface  interferes  with  some  other  phase  of  the  Reynolds  stress 
production  mechanism  (say,  for  example,  with  the  secondary  instability 
resulting  from  Che  inflectionary  profile);  this  would  allow  Che  eddies  Co 
grow  larger  before  bursting,  increasing  the  time  between  bursts,  resulting 
in  a  shift  of  Che  pressure  spectrum  toward  lower  frequencies. 

It  is  also  possible  Chat  some  as-yet-undeCermined  mechanism  is  respon¬ 
sible  for  the  initiation  of  Che  large  eddies  (and  hence  for  their  frequency 
of  occurence)  and  that  the  compliant  surface  interferes  with  this  in  some 
way.  Some  light  may  be  shed  on  this  question  by  experiments  on  Che  influ¬ 
ence  of  polymer  additives  on  Che  time  between  bursts. 

It  is  also  conceivable  that  Che  finite  thickness  of  the  compliant 
surface  is  significant,  since  the  time  for  transmission  and  reflection  of  a 
signal  through  Che  layer  is  of  the  order  of  the  growth  period;  the  reflec¬ 
ted  disturbance  could  interact  favorably  with  Che  eddy  at  a  later  phase  in 


Finally,  eddies  such  as  we  envision  also  arise  due  to  the  nonlinear 
interaction  of  two  oblique  progressive  waves.  It  may  be  that  such  progres¬ 
sive  waves  are  present  in  the  sublayer  (a  proper  orthogonal  decomposition 
theorem  applied  to  the  turbulent  field  gives  an  Orr-Somerfeld  equation, 
which  has  such  solutions;  see  Lumley  1967).  Such  waves  would  be  much  more 
strongly  affected  by  a  compliant  boundary. 


7.  References. 


Bushnell,  D.  M. ,  J.  N.  Hefner  and  R.  L.  Ash  (1977)  Effect  of  compliant  wall 
motion  on  turbulent  boundary  layers.  Physics  of  Fluids  Supplement  Vol.  20, 
pp.  S31-S48. 

Hansen,  R.  J.  ,  D.  L.  Hunston,  C.  C.  Ni,  M.  M.  Reischman  and  J.  W.  Hoyt 
(1980)  "Hydrodynamic  drag  and  surface  deformations  generated  by  liquid 
flows  over  flexible  surfaces”,  in  Viscous  Drag  Reduction  (G.  R.  Hough,  ed. ; 
Vought  Advanced  Technology  Center;  Dallas,  TX)  In  Press. 

Luroley,  J.  L.  (1964)  "The  reduction  of  skin  friction  drag",  in  Ship  Motion 
and  Drag  Reduction,  ACR-112,  pp.  915-946,  (Office  of  Naval  Research,  Dept 
of  the  Navy,  Washington,  D.  C.). 

Lumley,  J.  L.  (1967)  "Structure  of  inhomogeneous  turbulent  flows",  in 
Atmospheric  turbulence  and  radio  wave  propagation  (A.  M.  Yaglom,  V.  I. 
Tatarski,  eds.;  Moscow;  Nauka) .  pp.  166-188. 

Lumley,  J.  L.  (1971)  "Some  comments  on  the  energy  method",  in  Developments 
in  Mechanics.  Vol.  6,  pp.  63-88  (L.  H.  N.  Lee  and  A.  H.  Szewczyk,  eds.; 
Notre  Dame  Press,  Notre  Dame). 

Lumley,  J.  L.  (1977)  Drag  reduction  in  two  phase  and  polymer  flows.  Physics 
of  Fluids  Supplement,  Vol.  20,  pp.  S64-S71. 

Serrin,  J.  (1959)  On  the  stability  of  viscous  fluid  motions,  Archive  for 
Rational  Mechanics  and  Analysis,  Vol.  3,  pp.  1-13. 

Tennekes,  H.  and  J.  L.  Lumley  (1972)  A  First  Course  in  Turbulence  (MIT 
Press;  Cambridge,  MA). 

Willmarth,  W.  W.  and  T.  J.  Bogar  (1977)  Survey  and  new  measurements  of 
turbulent  structure  near  the  wall.  Physics  of  Fluids  Supplement,  Vol.  20, 
pp.  S9-S21. 


9 

o 

(T 

I 

I 

I 

I 


FLEXIBLE 


Figure  2.  Plot  .of  equation  5.15  for  a  fixed  value  of  normalized  growth 
rate  for  the  rigid  wall  case  (the  compliant  wall  case  is  virtually  indis- 
t inguishable) . 
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Figure  3.  Plot  of  critical  growth  rate  and  wavenumber  versus  Reynolds 
number  for  the  rigid  wall  and  compliant  wall  cases.  These  are  the  wavenum¬ 
bers  corresponding  to  the  minimum  value  of  Reynolds  numbers  for  each  growth 
rate  from  curves  such  as  that  shown  in  Figure  2.  Note  that  lambda  here  is 
as  defined  in  Lumley  (1971),  that  is,  the  growth  rate  normalized  by  sub¬ 
layer  variables. 
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Some  Comments  on  the  Energy  Method 
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ABSTRACT 

It  iw  suggested  that  the  large  orgamaed  motions  which  are  observed  in  inhomo¬ 
geneous  turbulent  flows  may  be  those  the  net  energy  of  which  as  the  largest 
growth  rate  (the  smaller  scale  turbulence  being  replaced  by  a  constitutive  rela¬ 
tion).  This  is  similar  to  the  classical  energy  method  for  stability  analysis;  a 
modern  extension  of  this  method  (by  Serrin,  1959)  is  described,  and  the  impli¬ 
cations  of  the  method  carefully  evaluated.  It  is  concluded  that  the  method  de¬ 
scribes  dynamically  possible  motions;  the  woric  of  Petrov  (1938),  in  w)uch  the 
contrary  conclusion  is  reached,  is  examined.  The  motions  are  found  to  be  nec¬ 
essarily  nonlinear,  and  to  evolve.  The  method  is  applied  to  longitudinal  rolls  in 
the  wall  region  of  boundary  layer  flow  (found  by  Joseph,  1966,  to  be  least  stable 
in  Couette  flow,  and  observed  in  the  wall  region).  The  resulting  equations  are 
solved  by  asymptotic  techniques;  these  make  clear  that  unstable  eddies  experi¬ 
ence  viscosity  only  near  the  wall;  that  streamwise  dist-  .bances  are  produced 
from  the  mean  gradient  by  vertical  disturbances;  that  in  the  viscous  region  tran¬ 
sverse  disturbances  decay,  while  still  producing  streamwise  disturbances.  As 
a  check  on  the  asymptotic  analysis,  it  is  also  applied  to  rectlinear  Couette  flow, 
where  it  gives  the  classical  value  for  the  zero- growth- rate  Reynolds  number 
(Serrin,  1959;  Joseph,  1966).  The  eigenvalue  relation  for  the  wall  layer  indi¬ 
cates  very  large  growth  rates  (in  sublayer  variables)  at  observed  Reynolds  num¬ 
bers,  in  agreement  with  observation.  Predicted  wavenumbers  are  larger  than 
observed;  a  way  is  suggested  by  which  the  reduction  in  shear  due  to  the  presence 
of  the  eddies  would  reduce  the  wavenumber  to  that  observed. 


INTRODUCTION 


It  is  now  well  known  that  inhomogeneous  turbulent  flows  exhibit  recurrent  struc¬ 
tures  in  the  velocity  field,  which  are  referred  to  as  "large  eddies"  (Townsend, 
1956).  The  dynamical  role  which  these  motions  play  was  discussed  at  length  in 
Townsend  (1956),  and  several  guesses  made  as  to  the  form  of  eddy  required  to 
explain  the  measurements  in  various  flows.  Suggestions  have  been  made  (L,um- 
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ley,  1967,  1970b)  as  to  how  these  motions  might  be  identified  objectively  in  cor¬ 
relation  data.  The  dynamical  equations  obeyed  by  these  eddies  have  been  dis¬ 
cussed  (Townsend,  1956;  Lumley,  1967).  Both  discussions  suggest  replacing  the 
turbulent  velocities  of  smaller  scale  by  a  constitutive  relation;  that  is,  treating 
the  large  eddy  deterministically,  as  weather  is  treated  meteorologically,  while 
the  smaller  scale  turbulence  is  treated  statistically.  Although  this  is  justified 
in  meteorology,  it  is  probably  not  in  Bows  in  the  laboratory:  it  is  not  hard  to 
show  that  a  spectral  gap  half  a  decade  wide  is  required,  for  such  a  treatment  to 
be  rigorously  justified  (Lumley,  1970c).  In  laboratory  Bows,  there  is  some¬ 
times  a  gap  of  sorts,  but  it  is  neither  wide  enough  nor  deep  enough.  Neverthe¬ 
less.  there  are  even  less  propitious  applications  of  this  technique  (to  cases  with 
no  gap  at  all)  which  have  been  remarkably  successful,  such  as  the  Heisenberg 
(1948)  treatment  of  the  spectrum.  Certainly  this  approach  does  not  violate  phy¬ 
sical  laws  (such  as  the  second  law  of  thermodynamics),  it  can  be  made  to  meet 
global  requirements  (such  as  dissipating  the  proper  amount  of  energy)  and  can 
be  endowed  with  certain  known  characteristics  of  the  real  situation  (such  as  vis¬ 
coelastic  behavior;  c.  f.  Lumley,  1970a);  hence  there  is  a  certain  justiBcation 
in  hoping  for  qualitatively  correct  results.  The  method,  in  addition  has  the  vir¬ 
tue  of  simplicity. 

The  equations  of  motion  for  the  large  eddy  are  now  reduced  to  the  Navier-Stokes 
equations,  if  a  simple  eddy  viscosity  is  assumed,  or  some  more  complicated 
equivalent,  if  a  more  realistic  constitutive  equation  is  adopted.  Although  this 
represents  a  considerable  simplification,  it  still  does  not  make  possible 
straight  forward  calculation:  the  number  of  known  exact  (or  even  asymptotic) 
solutions  to  the  Navier-Stokes  equations  is  remarkably  small.  Ordinarily,  for 
given  boundary  conditions,  the  variety  of  solutions  that  can  be  generated  by  var¬ 
ious  initial  conditions  is  considerable;  which  one  is  the  large  eddy?  According¬ 
ly,  another  principle  appears  to  be  needed,  to  isolate  the  large  eddy.  It  has 
been  suggested  (Lumley,  1967)  that  the  mean  velocity  profile  might  be  neutrally 
stable  (in  the  small  disturbance  sense)  to  the  large  eddy,  thus  identifying  the 
latter  as  the  eigenfunction  of  the  Orr-Sommerfeld  (or  equivalent)  equation  cor¬ 
responding  to  the  minimum  value  of  Reynolds  number.  Unfortunately,  direct 
calculation  (Reynolds,  Tiederman,  1967)  indicates  that  turbulent  boundary  layers, 
at  any  rate,  are  stable  to  all  such  disturbances. 

We  are  proposing  here  another  principle,  that  the  large  eddy  observed  is  that 
motion  which  can  most  efficiently  extract  energy  from  the  mean  motion  and  give 
up  as  little  as  possible  to  the  (turbulent)  dissipation.  This  is  the  motion  which 
can  therefore  grow  largest  under  given  circumstances.  Application  of  this 
principle  gives  a  well  defined  and  relatively  simple  eigenvalue  problem,  at  least 
for  the  case  of  a  Newtonian  constitutive  relation.  In  this  case,  it  is  formally 
equivalent  to  what  is  known  classically  as  the  energy  method  of  stability  analysis. 

We  will  restrict  ourselves  to  a  Newtonian  constitutive  relation,  since  we  have  in 
mind  applying  this  analysis  to  the  wall  region  of  a  boundary  layer  or  pipe  Bow. 

In  such  a  Bow,  the  maximum  shear  and  minimum  viscosity  is  in  the  sublayer, 
which  is  Newtonian;  one  may  reasonably  expect  that  the  dynamics  of  such  large 
eddies  as  exist  may  be  dominated  by  Newtonian  phenomena.  In  fact,  we  will 
Bnd  from  the  asymptotic  analysis  that  the  nature  of  the  effective  viscosity  out¬ 
side  the  sublayer  has  no  inBuence  on  the  energy  balance  of  the  eddies.  This 
would  certainly  not  be  generally  true;  in  attempting  to  calculate  the  form  of 
large  eddies  in  shear  layers,  for  example,  one  might  expect  the  character  of 
the  effective  constitutive  relation  to  be  important,  since  the  maximum  shear 
region  occurs  where  the  turbulence  is  strong;  the  energy  balance  then  will  be 
between  production  and  turbulent  dissipation,  and  the  effective  constitutive 
relation  is  certainly  not  Newtonian  (Lumley,  1970a). 

THE  ENERGY  METHOD  OF  SERRIN 


The  classical  energy  method  is  associated  with  Orr  (1907).  This  was  limited 
to  two-dimensional  disturbances  to  the  basic  Bow,  and  continuity  was  satisfied 
by  the  use  of  a  stream  function.  In  1959,  Serrin  extended  the  method  to  three- 
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dimensional  disturbancea,  using  Lagrange  multipliers.  Our  derivation  will  be 
restricted  to  a  Newtonian  constitutive  relation,  and  is  a  trivial  generalization  of 
Serrin's,  using  a  slightly  different  method. 

We  are  seeking  the  disturbance  u.  to  the  basic  flow  which,  subject  to  the  con¬ 
dition 

Ui_i  =  0  (1) 

has,  at  a  fixed  time  t,  for  a  fixed  value  and  distribution  of  v,  an  extremal  value 
of 

a  =  (l/2E)dE/dt  (2) 

where  E  is  the  global  disturbance  energy.  That  is,  the  disturbance  which  can 
most  efficiently  extract  energy  from  the  basic  flow,  losing  as  little  as  possible, 
will  have  the  largest  growth  rate.  That  this  extremism  is  truly  a  maximum  can 
be  seen  more  mathematically  by  considering  that,  for  a  fixed  distribution  of  v, 
the  dissipation(which  must  always  be  nonnegative)  may  be  made  arbitrarily  large 
by  selecting  u^  having  sufficiently  small  characteristic  length  scales.  The  pro¬ 
duction  is  bounded  (relative  to  the  energy)  above  and  below  (by  the  negative  of 
the  least  and  greatest  eigenvalues  of  the  strain  rate  tensor  of  the  basic  flow). 

Oue  to  the  dissipation,  however,  the  range  of  growth  rates  is  bounded  only  above, 
and  unbounded  below. 

Using  the  integral  of  the  equations  of  motion,  we  have  from  (2) 

a  =  [  -  JS.  .UjU.dV  -  jV2s.  .dVj  /;u.u.dV  (3) 

where  and  sjj  are  respectively  the  strain  rate  tensors  of  the  basic  flow  and 
of  the  disturbance.  Note  that  all  transport  terms  vanish,  since  they  are  conser¬ 
vative.  Applying  the  usual  techniques  of  the  calculus  of  variations .  w*  obtain 
(indicating  by  the  maximum  value) 


J{q'*’u.  +  S.jUj  .  2  (vs.j)  pSu.dV  =  0 


Now,  if  the  variation  is  over  incompressible  motions,  the  6u.  are  not  indepen¬ 
dent,  but  are  related  by  ^ 


S  u.  .  =  0 
1,  t 


We  may  specify  6U|  and  Su2  independently,  for  example,  and  then  6  uj  is 
given.  If  we  write 

o^Ui  +  S..U.  .  2  (vs.y  .  =  X.  (6) 

for  convenience,  then  (4)  becomes 

JX.Su.dV  =  J  (XjSUj  +  X2*“2  *  'I''’  =  “ 

If  we  set  X,  =  -  TT  3  /  p,  where  tt  is  an  arbitrary  function,  then  integration  by 
parts  (using^the  condition  that  u^  and  6u^  either  vanish  on  the  boundaries  or 
satisfy  a  cyclic  boundary  condition)  gives 

J'X.eu.dV  =  X,(Xj6Uj  +  X^SUj  +8uj  j  n/p)  dV 
=  J  (XjSuj  +  X^Su^  -  [6Uj  j  +  5'i2 


JfXj  +  TT_  j/p)«Uj+(X2+TT_2/p)8“2J 
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(9) 


Since  Suj  and  6u2  are  independent,  we  have 

®^“i^  Vj  '  .j  “i.i"  ° 

This  is  the  same  as  that  derived  in  Serrin  (1959)  with  the  exception  of  the  in¬ 
clusion  of  a  non- zero  growth  rate,  and  variable  viscosity.  The  pressure  n  is 
that  required  to  assure  (1),  and  hence  may  be  regarded  as  a  Lagrange  multiplier. 

This  has  been  applied  by  Serrin  (1959)  to  Couette  flow  with  uniform  viscosity  and 
a'*'  =  0.  Two-dimensional  disturbances  of  the  same  flow  were,  of  course,  cal¬ 
culated  by  Orr  (1907);  he  found  that  the  Reynolds  number  at  which  a'*'  =  0  (the 
critical  Reynolds  number)  for  such  disturbances  is  6.65,  based  on  the  shear 
velocity  and  the  half  wldt^  Although  Orr  (1907)  said  "Analogy  with  other  pro¬ 
blems  leads  us  to  assume  that  disturbances  in  two  dimensions  will  be  less  stable 
than  those  in  three;  ...  "as  pointed  out  by  Joseph  (1966),  Serrin  (1959)  found 
that  longitudinal  rolls  could  survive  to  a  lower  Reynolds  number,  namely  4.  54 
(on  the  same  basis);  evidently  such  disturbances  can  extract  energy  more  ef¬ 
ficiently.  Serrin  (1959)  did  not  show  that  this  type  of  disturbance  was  the  most 
efficient,  i.  e.  -  that  it  corresponded  to  the  largest  growth  rate  (=0)  at  a  given 
Reynolds  number,  or  to  the  lowest  Reynolds  number  at  a  given  growth  rate. 
Joseph  (1966)  however,  by  considering  the  combined  problem  of  Couette  flow  of 
a  Boussinesq  fluid  heated  from  below,  showed  that,  in  fact,  the  longitudinal  rolls 
are  the  most  efficient. 

Taylor  (1960)  suggested  using  this  type  of  analysis  for  a  purpose  different  from 
ours;  namely,  viewing  a  Couette  flow  as  two  viscous  sublayers  face  to  face,  the 
thickness  Reynolds  number  (defined  as  above)  corresponding  to  a  maximum 
growth  rate  of  zero  appears  to  give  the  thickness  of  the  layer  next  to  the  wall  in 
which  no  disturbance  can  survive  without  importing  energy  from  above.  He  used 
the  analysis  of  Lorentz  (1907)  who  did  not  actually  solve  the  extremum  problem, 
but  guessed  at  a  solution.  Using  Serrin  and  Joseph's  solution,  tliis  would  sug¬ 
gest  a  viscous  sublayer  thickness  of  y'*'  *  4.  54,  wliich  is  near  the  point  at  wliich 
the  mean  velocity  profile  bends  away  from  the  linear,  although  well  below  the 
point  (~9)  at  which  measured  dissipation  first  exceeds  production  (Townsend, 
1956), 

Before  applying  (9)  to  a  specific  calculation,  we  must  examine  the  implications 
of  the  method,  so  that  we  will  know  what  to  expect  of  it, 

CRITIQUE  OF  THE  ENERGY  METHOD 

The  energy  method,  as  applied  to  flows  without  buoyancy  ,  has  been  subject  to 
extensive  criticism.  For  example,  Lin  (1955)  states:  "Even  here  only  the  lower 
limit  for  the  critical  Reynolds  number  can  be  expected,  because  stability  must 
be  estabiished  for  all  disturbances  in  the  present  method,  while  in  reality  only 
those  satisfying  the  hydrodynamic  equations  need  to  be  considered.  The  inclu¬ 
sion  of  the  spurious  disturbances  ..."  (pag  :  59)  and  "The  energy  equation  thus 
gives  critical  Reynolds  numbers  that  are  too  low,"  (page  61).  Serrin  (1959) 
states:  "It  is  important  to  note  that  the  energy  method  cannot  provide  accurate 
knowledge  of  the  limits  of  stability,  such  as  can  be  gained  from  the  linearized 
perturbation  theory  ...  The  reason  is  that  in  the  energy  method  one  establishes 
stability  relative  to  arbitrary  disturbances,  while  in  reality  only  those  satisfy¬ 
ing  the  hydrodynamical  equations  need  be  considered,  "(page  4).  Joseph  (1966) 
states  "The  functions  which  are  admitted  for  review  need  not  be  possible  solu¬ 
tions  of  the  conservation  equations  (momentum  and  energy).  It  is  in  this  sense 
that  one  may  speak  of  dynamically  inadmissible  disturbances,"  (page  181).  And 
further:  "Evidently  ...  the  energy  method  reflects  a  sen  itivity  to  spurious  and 
dynamically  inadmissible  disturbances."  (page  182).  Non.>  of  these  allegations 
is  supported.  However,  Monin  and  Yaglom  (1971)  state  "Analyzing  this  fact 
[that  the  energy  method  gives  lower  critical  Reynolds  numbers  than  the  linear 
theory]  Petrov  (1938)  came  to  the  conclusion  that  the  values  of  t  at  which 

*with  e  0,  and  uniform  viscosity. 
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F  ['{r]  takes  a  maximum,  taking  into  account  the  time  variation  of  all  the  func¬ 
tions,  will  apparently  in  no  case  generate  a  dynamically  possible  motion.  Thus 
the  energy  method  can  never  give  an  exact  value  of  Re^^  min’  only  suit¬ 

able  for  making  preliminary,  very  rough  estimates  of  tUs  value.  " 

The  question  of  why  the  critical  Reynolds  numbers  predicted  by  the  energy  met¬ 
hod  in  flows  without  buoyancy  are  substantially  lower  than  those  predicted  by  the 
linear  theory'’’  is  certainly  an  interesting  one,  although  whether  they  are  "too 
low"  as  stated  by  Lin  (supra)  and  in  what  sense,  remains  to  be  seen.  The  ex¬ 
planations  given  by  Lin,  Serrin  and  Joseph  are  clearly  not  satisfactory,  since 
functions  are  "adrnitted  for  review"  at  a  fixed  time  only;  any  velocity  field  at  a 
fixed  time  satisfies  the  hydrodynamical  equations,  which  serve  to  determine  its 
evolution,  i.  e.  -  the  time  derivative.  The  only  substantial  argument  appears  to 
be  that  of  Petrov  (1938). 


Since  the  paper  of  Petrov  (1938)  appears  never  to  have  been  translated,  a  trans¬ 
lation  of  the  relevant  section  (pages  20-24)  is  appended  (some  obvious  misprints 
in  the  equations  have  been  corrected).  The  remainder  of  the  Petrov  paper  is  a 
concise  introduction  to  classical  small  disturbance  stability  theory  and  the  ener¬ 
gy  method,  and  need  not  concern  us  here.  Petrov,  of  course,  is  reasoning  pri¬ 
marily  about  linearized  disturbances  having  exponential  behavior  in  time. 

His  argument  may  be  summarized  as  follows:  the  small  disturbance  equations 
obtained  from  the  Navier-Stokes  equations,  and  the  equations  obtained  from  the 
energy  method  are  not  the  same;  both  cannot  be  satisfied  simultaneously  (the 
difference  is  Petrov's  equation  (22));  hence,  disturbances  considered  in  the 
energy  method  do  not  satisfy  the  (small  disturbance  limit  of  the)  Navier-Stokes 
equations. 

This  conclusion  is  correct  when  applied  to  the  small  disturbance  limit  of  the 
Navier-Stokes  equations  (though  misleading);  it  is  not  correct  for  the  full  equa¬ 
tions.  The  difference  lies  in  the  interpretation  of  the  time  derivative.  The 
energy  method  specifies  a  global  growth  rate;  locally  the  growth  rate  may  be 
greater  or  less,  with  the  non-linear  terms  transferring  energy  from  one  region 
to  another  as  required.  These  non-linear  transport  terms,  being  conservative, 
do  not  appear  in  the  integrals,  so  that  the  equations  of  the  energy  method  apply 
to  large  ^sturbances  whether  this  is  intended  or  not.  The  difference  between 
the  equations  of  the  energy  method  and  the  Navier-Stokes  equations  serves  to 
determine  the  local  growth  rate.  Only  in  the  linearized  case,  where  no  mechan¬ 
ism  is  present  to  transport  energy  from  one  level  to  another  (in  a  parallel  Qow) 
is  the  growth  rate  required  to  ce  the  same  at  all  levels,  and  only  then  does  in¬ 
consistency  result.  Hence,  it  is  the  assumption  of  linearity  that  is  inconsistent. 
Let  us  examine  the  equations. 


The  condition  (2)  does  not  imply  that  8u./3t  =  au-  at  each  point  in  the  fluid  at  the 
instant  in  question.  The  growth  rate  of  the  disturbance  may  be  larger  or  smaller 
than  optimum  locally,  so  long  as  (2)  is  satisfied.  In  fact,  of  course,  must 
also  satisfy  the  Navier-Stokes  equation 


u.  +  U.  .u.+u.  ,U.+  u.  .u  .  =  -  p  ./p  +  (v(u.  .  +  u.  .)  )  . 

1  i.J  }  i.J  J  i.J  }  ^,1  '  i.j  j,i'  ',j 

and  the  local  growth  rate  can  be  obtained  by  substracting  (9)  from  (10) 

u.-a*u.  =  -  n.  .u.-u.  .U.-u.  .u.  -  (p-tt)  .  /p 
1  I  “ij  j  I.J  J  I.J  J  '.1 


(10) 


(11) 


In  Joseph  (1966)  it  is  shown  that  in  several  flows  with  buoyancy,  the  critical 
Reynolds  numbers  predicted  by  the  energy  method  are  close  to,  or  the  same  as, 
those  predicted  by  the  linear  theory. 
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where  U-  j  =  S^j  +  Cly,  being  the  antisymmetric  part.  As  required  by  (2), 
multiplication  dl  (llTby  ^  and  integration  over  the  region  gives 

dE/dt  =  (d/dt)  j’u.UjdV/a  =  Ju.ii.  dV  =  o+Ju.u.dV  (12) 

If  our  explanation  is  correct,  equation  (11),  when  restricted  to  a  linearized  two- 
dimensional  disturbance  of  a  two-dimensional  parallel  flow,  should  give  Petrov's 
equation  (22).  Let  U.  =  (U(X2)>  0,  0). 

0  u'/2  0  0  U /2  0 

S..  =  u'/2  0  0  ,  a.  =  -U  /2  0  0  (13) 

0  0  0  0  0  0 
I 

where  U  =  U  ^  satisfy  (1)  by  the  use  of  a  stream  function  Uj  =  <l',2’“2 

X-  j.  EUminating  the  pressure  between  1  and  2  components  of  the  linearized 
(11)  gives 

(a/at)7^»  -  =  u"  f _  j/2  -  U7^  ♦,  1  *  u'  21 

Now,  in  the  linear  case  (considering  constant,  uniform  viscosity)  the  coefficients 
of  (10),  or  of  the  equation  for  the  stream  function  obtained  from  it,  are  indepen¬ 
dent  of  ,  t,  and  hence  the  solution  may  be  written  as  a  sum  of  terms  of  the 
form 

♦  =  ^  e‘®<*l-‘=‘)  c  =  c^  +  ic.  (15) 

where,  from  (2)  and  (15)  we  necessarily  have 

a*  -  acj  (16) 

Substituting  this  in  (14)  we  obtain  immediately  Petrov's  equation  (22)  • 

(U  -  c^)  («"  .  ah)  -  v'h/Z  +  u'  ®'  =  0  (17) 

That  this  is  inconsistent  is  clear  because,  as  Petrov  points  out,  one  parameter 
(the  viscosity)  is  absent  from  (17);  hence  an  eigen-solution  of  (17)  cannot  be  con¬ 
sistent  with  an  eigen- solution  of  the  Orr-Sommerfeld  equation  (obtained  from  the 
linearized  equation  for  the  stream  function  obtained  from  (10),  on  substitution 
of  (15)  ),  since  the  latter  contains  the  viscosity. 

Hence,  we  must  conclude  that  the  extremal  disturbance  generated  by  the  energy 
method  is  necessarily  non-linear,  entailing  transport  of  disturbance  energy, 
and  necessarily  has  a  local  growth  rate  of  the  disturbance  energy  differing  from 
the  global  growth  rate,  but  averaging  spacially  to  the  latter.  The  disturbance 
does  satisfy  the  hydrodynamical  equations,  which  serve  to  determine  the  local 
growth  rate.  Since  the  class  of  disturbances  considered  is  no  smaller  than,  and 
may  be  larger  than  that  considered  in  the  linear  theory  (which  necessarily  have 
no  transport,  and  the  same  growth  rate  everywhere)  the  Reynolds  number  ob¬ 
tained  for  the  extremal  value  o'*'  =  0  is  no  larger  than,  and  may  be  less  than, 
that  obtained  from  the  linear  theory. 

As  an  aside,  it  may  be  noted  that,  since  the  (necessarily)  non-linear  disturbance 
generated  by  the  energy  method  for  a  =  a*  (necessarily)  has  local  growth  rates 
different  from  o^  it  will  evolve  in  form;  at  an  instant  following  (or  preceding) 
the  instant  of  analysis,  it  will  have  changed  and  may  no  longer  be  optimum. 

But,  if  the  optimum  disturbance  evolved  continuously  through  optimum  states, 
equation  (9)  would  have  to  be  satisfied  in  some  neighborhood.  We  would  thus 
have  two  eigenvalue  problems  in  this  neighborhood:  either  (9)  and  (10),  or 
equivalently  (9)  and  (11).  But  we  may  use  Petrov's  argument  that  the  sets  (9) 
and  (11)  cannot  be  simultaneously  satisfied,  due  to  the  absence  of  the  viscosity 


in  (11).  Hence,  (9)  end  (11)  can  at  most  be  simultaneously  valid  at  a  point,  and 
the  optimum  disturbance  must  evolve  to  one  w)uch  is  not  optimum,  and  hence  has 
a  growth  rate  less  than  the  extremal  value  a*.  Hence,  the  disturbance  energy 
has,  at  the  instant  of  analysis,  a  global  growth  rate  of  Za* ,  and  in  some  neigh¬ 
borhood  of  that  time  a  smaller  growth  rate.  At  no  time  can  the  energy  have  a 
larger  growth  rate  than  the  optimal  one;  if  the  optimal  disturbance  were  not 
unique,  the  disturbance  energy  might  at  later  times  evolve  to  other  optimum 
states,  producing  points  of  growth  rate  equal  to  2a'*';  by  the  preceding  argument, 
however,  these  must  be  isolated,  so  that  the  disturbance  energy  has  a  growth 
rate  almost  everywhere  less  than  Za* .  Hence,  the  optimal  disturbance  energy 
is  bounded  from  above  by  the  exponential  with  growth  rate  2a'*'  (  a  result  of 
Serrin,  19S9,  obtained  in  a  different  way),  and  has  this  growth  rate  at  the  in¬ 
stant  of  analysis;  at  other  times,  the  growth  rate  is  almost  everywhere  smaller. 
Thus,  optimum  disturbances  obtained  for  a^  =  0  have  a  disturbance  energy 
which  is  monotone  decreasing  almost  everywhere. 

We  must  note  that  we  have  not  excluded  the  possibility  that  disturbances  which 
have,  at  the  instant  of  analysis,  positive  global  grow^  rates,  will  ultimately 
decay.  These  would  also  be  counted  stable.  For  tliis  reason  also  (and  perhaps 
primarily)  the  critical  Reynolds  number  given  by  the  energy  method  for  a**^  =  0 
is  conservative:  it  assures  that  no  disturbance  can  grow  ultimately,  and  that  at 
a  higher  value  some  disturbance  can  grow  initially;  it  does  not  guarantee  that 
there  is  a  disturbance  which  can  grow  ultimately  at  a  higher  Reynolds  number. 
Thus,  one  says  that  at  a  given  Reynolds  number  the  energy  method  determines 
stability,  but  cannot  determine  instability.  In  fact.  Joseph  has  recently  proved 
(1971)  that  longitudinal  rolls  in  a  shear  flow,  although  they  have  an  imtially  in¬ 
creasing  disturbance  energy,  always  decay  ultimately. 

We  may  speculate  that,  in  flows  with  heat  transfer,  there  is  a  smaller  disparity 
between  the  linear  and  non-linear  critical  values  because  the  buoyancy  provides 
a  mechanism  for  transport  of  disturbance  energy  normal  to  the  flow  even  in  the 
linearised  case. 

APPLICATION  OF  THE  METHOD  TO  THE  CONSTANT  STRESS  LAYER; 

THE  ASYMPTOTIC  SOLUTION 


We  are  going  to  apply  this  method  to  the  constant  stress  wall  layer  of  a  turbulent 
boundary  layer  or  channel  flow  with  sero  pressure  gradient.  Rather  than  at¬ 
tempt  to  find  the  general  type  of  disturbance  which  is  most  efficient,  we  will  use 
the  result  of  Serrin  (1959)  and  Joseph  (1966)  for  the  Couette  flow,  presuming 
that  here  also  longitudinal  rolls  will  be  most  efficient;  thus  we  will  assume  a 
disturbance  independent  of  xj  (the  streamwise  direction)  and  periodic  cross¬ 
stream.  Setting  Uj  *  (U(x2),  0,  0). 


ikx, 

Uj  =  u  e  3 


•I-,,  ikx,,  ikx- 

Uj  =  ik(r  e  3,  Uj  =  -  I)  e  3 


(18) 


where  u,  and  i  are  functions  of  x^.  and  a  prime  denotes  differentiation  with 
respect  to  Xj,  equation  (9)  becomes^ 

c^u  +  U  iki|;/2  =  (vu  )  -  It^  vu 

ika^*  +  U  u/2=-n/p+2  (vik*  )  -  ikv(i)  +  k^t)  (19) 

-  a* i  =  -  TTik/p  -  {v(t  +  k^it)]  +  2vk^tr 
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Eliminatine  the  pressure  tt  between  the  second  and  third  of  (19)  results  in 
(writing  sD^-k^,  D  =  d/dy) 


+  U  iku/2  =  (v7^7^ 


+  2v  7  D  +  V 


[  +  k^]  )  ♦ 


O^u  T  U  iki|t/2  =  (v7^ 


+  V  D)  u 


(20) 


In  the  coni^tant  stress  layer  of  a  turbulent  boundary  layer  we  may  write 

N>  3  u<r  ^/U  ,  where  u^  is  the  shear  velocity.  Thus,  in  the  viscous  sublayer, 


V  «  V 


m 


the  molectilar  value,  while  in  the  logarithmic  region  we  have 


(21) 


V  ~  u^Ky 


(22) 


where  K  is  von  Karman's  constant,  ~  0.4.  We  may  use  Reichardt's  (1951) 
expression,  which  agrees  with  measurement  reasonable  well  through  the  entire 
region  (see  Figure  1): 

v  =  Vj^[l  +  iC  (y'*'  -  R  tanh  y  ■*' /R)I  (23) 

where  R  is  the  point  where  the  extended  wall  profile  meets  the  extended  logar¬ 
ithmic  profile,  roughly  11.6  and  y'*'  =  yu^/v^.  We  will  find  that  our  final  ans¬ 
wer  is  only  weakly  dependent  on  the  particular  form  chosen  for  (23);  it  is  only 
essential  A^t  the  first  two  derivatives  vanish  at  the  wall.  Using  (23),  it  can  be 
seen  that  v  is  zero  in  the  sui^^ayer,  rises  near  R,  and  takes  the  constant  value 
K  in  the  logarithmic  region;  v  is  zero  in  the  sublayer  and  in  the  logarithmic 
region,  be^ng  non-zero  (and  positive)  only  near  y'*'  =  R.  From  the  relationship 
V  3  u.^  evidently  express  the  mean  velocity  gradient  also  in  terms 

of  the  function  (23).  If  we  non-<fimensionalize  y  by  the  value  corresponding  to 
y*  =  R-  say  y*  so  that  the  dimensionless  y  =  y+  /R  s  y/y*,  and  with 
\  3  y,^  ^  o'*’  /vnj,  then  the  equations  can  be  written  as 


X7^t  +  (R^k/2f)  u  3  17^7%  +  2f'  7^D*  +  f"  (D^  +  k^)  * 

2  2  ' 

Xu  +  (R  ik/2f)  f  3  f7  u  +  f  Du 

with  Reichardt' s  form  for  f 

f  3  1  +K  R(y  -  tanh  y  )  (25) 

We  see  (c.  f.  Figure  2)  that  the  term  in  R^  is,  for  R  =  11.  6  (the  true  value  in 
the  boundary  layer)  nearly  10^  times  the  terms  on  the  right-hand  side  near 
y  3  1.  As  y  increases,  ^e  terms  do  not  become  comparable  until  y  ~  10^. 
This  suggests  an  asymptotic  solution  (Cole,  1968)  using  R^  as  the  large  para¬ 
meter  (but  holding  the  implicit  R  which  appears  in  the  definition  of  f  fixed*). 
In  order  to  make  a  sensible  (non-trivial)  problem,  we  must  assume  that  X  = 
0(R‘)  also;  that  is,  that  outside  the  (viscous)  wall  layer  production  and  growth 
balance,  while  in  the  wall  layer,  viscous  loss  is  important. 


*It  appears  to  be  possible  to  construct  an  expansion  allowing  this  R  to  vary 
also;  it  is  more  complicated,  however,  in  that  three  distinguished  limits  are 
found,  and  the  behavior  of  the  solutions  in  our  range  of  parameters  is  the  same. 
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Suificiently  far  outside  (for  y  >  10^)  there  must  again  be  a  region  where  viscos¬ 
ity  is  important,  but  there  the  solution  i|  already  so  small  that  we  may  safely 
ignore  the  transition  to  this  new  region*  .  Thus,  writing  g^=X/R^,R=  e'^,  we 
have 


0^7^*  +  (ik/2f)  u  =  (fV^V^'l  +  Zf'y^Dm  +  f"(D^  +  k^)  i|,) 


8^  u  +  (ik/2f)  ii  *  e^(f7^u  +  f  u  ) 


(26) 


Now,  for  the  inner  problem  (next  to  the  wall)  we  can  write  n  =  y/e,  and 


I  =  Vq  (e)  (gjj  (h)  +  e^gj(n)  +  •••  ) 
u  =  Vq  (e)  (hjj  (t))  +  (n)  +  •  •  • ) 


(27) 


The  order  of  the  boundary  layer  thickness  was  determined  from  the  necessity 
of  keeping  the  highest  order  derivative  in  the  equations.  The  order  of  the  sec¬ 
ond  order  terms  in  (27)  was  determined  from  the  desire  to  obtain  a  gj,  hj 
different  from  gn>  hg.  If  necessary  for  the  matching  (it  will  not  be)  another  gQ 
and  hg  of  order  between  and  1  might  be  added  in  (27).  The  equality  of  the 
orders  of  the  leading  terms  arises  from  the  need  to  retain  (i  terms  in  the  u 
equation,  and  u  terms  in  the  ^  equation.  The  equations  are  (after  cross- sub¬ 
stituting) 


nil  2 
«0  *  ®  «0 


=  0;  g 


S^8i 


=  ikhg/2  +  2k  g^ 


(28) 

l»0  -  ikgjj/2;  hj  -  g^hj  =  k^hg  +  ikgj/2 

I 

The  boundary  conditions  at  the  wall,  of  course,  are  1(0)  =  i|i  (0)  =  0;  u(0)  =  0. 
Applying  these,  the  solutions  which  we  need  become 

gg  =  D  (1  -  gn  -  e 

hg  .  D  [e*^’^  .  1  +  gn+  (gT,/2)e'®’^Jik/20^ 


gjaB(e'®"-l)+Bn[B+(<?/4+48V^D/24B^]-(l-4B^  (B^-n^-B^T!^/3>k-^D/8B^ 


-Bti  e'®^(9/4  +  4B^k^D/24B^  +  g^n^  e*®’V^D/320^  (29) 

We  have  also  discarded  terms  that  are  transcendentally  large  as  r  -*  >,  whicn 
would  preclude  matching.  The  constant  B  must  be  found  in  terms  of  D  by  the 
matching.  O,  of  course,  must  remain  arbitrary,  since  the  equations  are  homo¬ 
geneous.  In  addition,  ot  course,  we  will  have  to  determine  vg  from  the 
matching. 

It  is  interesting  to  note  in  equation  (28)  that  to  zeroth  order  (i  (i.  e.  -  the  lateral 
motioh)  simply  experiehces  a  balance  between  decay  and  dissipation;  the  stream- 
wise  motion,  however,  experiences  some  production  due  to  the  lateral  motion. 

1 

That  is,  we  will  find  that  it  is  the  matching  to  the  solution  between  1  and  10 
.that  determihes  the  eigenvalue  relations  to  the  order  in  which  we  are  interested. 
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Thia  agrees  with  the  conclusions  ol  Joseph  (1971).  From  tne  equation  for  hg  it 
IS  Clear  that  tne  vortices  are  sweeping  slow  moving  fluid  up  from  tne  wall  (and 
vice  versa)  producing  disturbances  in  the  streamwiae  motion. 

We  may  now  attack  the  outer  problem.  Since  the  an>plitude  of  the  solution  is  ar¬ 
bitrary,  we  may  take  it  to  be  of  order  unity,  writing 


♦  =  Gg  (y)  +  c^Gj  (y)  +  ••  • 
u  *  Hq  (y)  +  (y)  +  ••• 


(30) 


where,  again,  the  order  of  the  second  term  is  fixed  by  requiring  that  G,  be 
different  from  Gg  (a  term  of  intermediate  order  similar  to  Gg  could  be  in¬ 
cluded,  but  will  not  be  necessary).  To  first  order  we  have 


+  (ik/2f)  Gq  =  0 
S^7^Gq  +  (ik/2f)  =  0 


(31) 


The  first  of  these  is  the  "mixing  length"  assumption,  frequently  made  (see  Bake- 
well,  Lumley,  1967):  that  the  perturbations  in  u  are  those  produced  from  the 
mean  velocity  gradient  by  the  lateral  motion.  Combining  the  two  equations  (31), 
we  obtain 


Gq  -  k^  (1  -  l/4fV)  Cq  =  0  (32) 

It  is  fairly  simple  to  show  that  this  solution  cannot  match  with  (29)  unless  there 
is  a  turning  point  •  i,  e.  -  unless  1  -  l/4f^B^  has  a  zero  somewhere.  This  is 
a  simple  kinematic  requirement;  if  (32)  has  no  turning  point,  the  lateral  vel¬ 
ocity  (outside  the  viscous  layer)  is  monotone  with  distance  from  the  surface, 
producing  upward  and  downward  streaming  jets,  rather  than  cells.  These  do 
not  satisfy  the  requirement  that  the  stream  function  and  its  derivative  must  van¬ 
ish  at  infinity.  If  (32)  has  one  turning  point,  we  have  one  cell  (see  Figure  3). 

Two  turning  points  correspond  to  two  cells,  each  on  top  of  the  other,  and  so  on. 
Clearly,  more  cells  means  a  smaller  length  scale,  and  hence  greater  dissipation 
Thus,  the  cases  of  two  and  more  turning  points  correspond  to  a  lower  growth 
rate.  We  will  consider  only  the  case  of  a  single  turning  point. 

A  simple  form  of  the  solution  for  (32)  may  be  obtained  (Cole,  1968)  formally 
valid  for  large  k  (but  probably  satisfactory  for  k  as  small  as  3)  by  the  two  var¬ 
iable  expansion  procedure.  If  we  designate  by  y^  the  point  where 

28^f(y„)  =  1  (33) 

c 

then  for  y  >  y^.  we  can  write 

G.  =  [1- l/4f^B*J  +,  [1  -  l/4f^8*]^^^dy'  (34) 

Y  c 

where  we  have  discarded  the  positive  exponential.  Now,  for  y  <  y,.  the  solu¬ 
tion  takes  the  form  of  trigonometric  functions.  The  proper  branch  (i.  e.  -  the 
continuation  of  (34)  may  be  selected  by  constructing  an  inner  expansion  of  (32) 
near  the  turning  point,  as  is  discussed  in  Cole  (1968)  (see  Figure  4).  The  re¬ 
sult  is  for  y  <  y_, 
c 


72 


Gq  =  2  /2/  3  (sin  n/ 124- cos  tt/  12)  [1/ .  1  ]  ”^^^sin{3TT/4-kY  ) 

(35) 

=  yl  [1/48V  -  1]  dy' 

For  smaller  k,  the  coefficient  in  (12)  becomes  a  series  in  k'^,  the  terms  involv¬ 
ing  integrals  and  derivatives  of  1/48 'f^  -  1.  A  little  manipulation  of  the  equa¬ 
tions  is  sufficient  to  show,  however,  that,  as  written  in  (35),  the  coefficient  has 
vanishing  first  and  second  derivatives  at  the  wall,  while  the  full  series  repre¬ 
sentation  of  the  coefficient  has  vanishing  first  through  third  derivatives,  due  to 
the  vanishing  of  the  first  two  derivatives  of  f.  (Substitution  of  a  form  such  as 
(35)  with  arbitrary  coefficient  into  (32)  gives  an  equation  for  the  coefficient, 
from  which  the  derivatives  at  the  wall  may  be  determined;  if  the  first  derivative 
vanishes,  the  first  three  derivatives  vanish.  But  a  non- vanishing  first  deriva¬ 
tive  produces  a  term  which  could  only  be  matched  with  (29)  if  f'  or  f  had  been 
assumed  not  to  vanish  at  the  wall. )  Thus,  although  the  proper  form  of  the  coef¬ 
ficient  for  finite  k  will  change  the  shape  of  Cq  away  from  the  wall  somewhat, 
it  will  not  affect  the  matching  through  third  order;  we  will  need  only  second 
order. 

The  second  order  equations  for  and  Hj  are,  of  course. 


8^7^J  +  ikHj/2f  =  fy^V^Gg  +  2f'7^DGg  +  f"  (D^  +  k^)  Gg 
8^Hj  +  ikGj/2f  =  f7^Hg  +  fHg' 


(36) 


t 

and  the  right  handsides  may  be  evaluated  in  terms  of  Gq  and  Gg  using  equa¬ 
tions  (31)  and  (32).  We  are  primarily  interested  in  the  behavior  of  the  solutions 
near  the  wall,  however.  If  account  is  taken  of  the  fact  that  Gg(0)=0(e) ,  Gg  (0)  = 
0(1)  (which  we  can  only  know  from  the  first  order  matching;  i,  e.  that  the  invis- 
cid  boundary  condition  is  nearly  satisfied,  while  the  no  slip  condition  is  not)  then 
the  leading  term  (through  linear  terms  in  y)  is 


B^7^Gj  +  ikHj/2  a  (k^/4B'‘)^Gg 

8^Hj  +  ikGj/2  a  (ik^/SS^Gg 


(37) 


or 

g”  +  k^  (1/48^  -  1)  Gj  =  (2kVl68*°)  Gg  (38) 

To  linear  terms,  then,  the  particular  solution  will  be 

(Gj)p  »  (2k^/168^°)  (1/48 -  1)  ’‘Cg  (39) 

This  result  is  the  mathematical  expression  of  the  physical  trutn  that  near  tne 
wail,  the  variation  of  viscosity  with  y  is  not  important. 

The  general  solution  will  be  exactly  like  (35)  -  we  can  thus  write  to  tne  order  of 
linear  terms  in  y, 

i|r  *  Gg  +  e^[C  Gg  +  (2k^/16e‘°)  (1/48*  -  D'^Gg)  +  •••  (40) 
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where  C  is  an  unknown  constant  that  must  be  determined  from  the  matching. 

We  may  now  proceed  to  the  matching,  using  an  intermediate  limit  in  which  y/6 
“yj  =  0(1),  6(e)  0,6  /  e-^  •.  so  that  t)  =8yj  /e  The  inner  solution  be¬ 

comes  (we  will  do  the  matching  only  in  i) ) 

(^QD(l-86yj2e)+e^v-C-B+g(6yj/e)  (B+(9/4+4e'*)k^D/24e^] 

(41) 

-(l-4@'*)  [e^(6yj/€)^-e^  /3lk^D/8e^]  +•  •  • 

neglecting  transcendentally  small  terms. 

Writing 

3n/4  -  kY_  =  n(e)  +  ^1/46“*  -  1)  +  0  (6^  (42) 

near  the  wall,  we  can  write  (to  third  order) 


Gj,  =  P  sin  [u(e)  +  k(l/4e‘‘  -  1)  ^/^6yj  +  •••  ]  (43) 

where  P  is  a  constant.  Thus,  through  terms  0(6^),  0(^^), 

Go=P{M'U^/6+6yjk(l/4e^-l)‘/^(l.M^/2).k^(l/4p^-l)(u.U^/6)(6yj)^/2 

(44) 

.  k^  (1/46^- 1)^/^  (l-a^/3)  (6yj)V6  +  •••  } 
and,  from  (40) ,  i)  is  given  by 


^  =  (44)+Pe^CC+(2k^/16g^°)  (1/4B^-1)‘‘} 
{n-u^/6+«yjMl/49'*-l)‘'^  (l-uVz)  +  •••  } 


(45) 


where  the  second  term  is  correct  to  0(6). 

Matching  the  term  of  0  (6)  in  (45)  we  find  that  we  must  have  Vg  =  e  , 

-go  =  Pk(l/4g^-l)*^^  (46) 


Writing 


2  3 

u(e)  =  «Uj«  +  uije  +  uijS  +••• 


we  may  match  the  term  of  0  (e)  in  (41),  giving 

D  =  Pis, 


(47) 


(48) 


Since  there  is  no  term  of  0  (e^)  in  (41),  is?  =  0,  Matching  terms  of  O(e^),  we 
have 

-B  =Pl«)3  -  P/6  +  P(C  +  (2k^/16g*°)  (l/4g^-l)'‘]  lUj  (49) 
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Terms  of  0(e^6)  give 

@  [B  +  (9/4+  4e^)k^D/24e^I  =  -  Pk(l/4e‘*  - 

+  PtC+  (2k^/16e^°)  (l/49‘*-l)‘^]k(l/4e‘*-l)‘''^ 
Terms  of  0  (e6  give 


(50) 


-(1  -  49^  9^k^D/89^  =  -  Pk^{l/49‘‘-l)ii,j/2  (51) 

and  of  0  (6^) 

+  (l-49'*)9^k^D/24e^  =  -  Pk^(l/49‘*-l)^‘'^/6  (52) 

Aa  a  result  of  (46),  (51)  and  (52)  are  satisfied  exactly.  ^Equations  (49)  and  (50) 
represent  two  equations  in  three  unknowns;  fortunately  two  of  the  unknowns  occur 
in  the  same  combination,  giving 


uij  =  -  (1/48'*  -  1)  (1/489*+  l)k^/28^  (53) 

The  other  combination  of  constants  need  not  concern  us. 

We  can  thus  write,  using  (46),  (47),  (48)  and  (53),  our  eigenvalue  conditions: 

3TT/4.kY_g=-(l/49*.l)‘''^ek/9-(l/4e*-l)‘'^^(l/48e*+l)e^k^/28^+  •••  |54) 

to  0(e®),  where  Y  ^  is  given  by 

^-0  =  o/'  '‘y  '55) 

Expression  (54)  gives  the  offset  (see  Figure  5)  by  which  the  outer  (inviscid)  so¬ 
lution  fails  to  meet  the  inviscid  boundary  condition  to  leave  room  0(c)  for  the 
viscous  (eddy)  boundary  layer.  Equation  (54)  has,  in  general,  two  positive  roots 
for  k  (for  fixed  e,  9)!  vre  are  interested  in  the  critical  condition  when  there  is 
only  one  (double)  root.  This  may  be  obtained  by  differentiating  the  expression 
(54)  with  respect  to  k,  requiring  that  the  derivative  vanish  also  at  the  same  value 
of  k.  If  we  indicate  by  k+  tliis  critical  value,  we  obtain  (noting  that  S/c  =  /\  i 


(jnr l\)^  =  (4/3n)  (l/4B*-l)*^^(l/486*+l) 

(56) 

Y  A“  =  (1/48*-!)*''^ +(l/48*-l)^^5,i/48g4^j^l/3j3^2)/5^(j^/g, 2/3 

"0 

U  now  a  value  of  g  is  selected,  Y.q  may  be  calculated,  as  well  as  and 

/X  Y.q;  from  the  latter  then  we  have  /\  ,  and  from  the  former.  k+.  while 
from  R  =  we  obtain  R. 


If  Figure  6  we  show  a  plot  of  X  and  k^  versus  R  (holding  the  R  in  f  fixed  at 
11.  6).  It  is  worth  noting  that  the  form  of  the  viscosity  outside  the  sublayer  did 
not  influence  the  solution;  the  only  property  we  used  was  the  vanishing  of  the 
first  two  derivatives  at  the  wall,  a  kinematic  requirement  in  a  zero  pressure 
gradient.  The  form  of  the  mean  shear  only  influences  expression  (55). 
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COUETTE  FLOW 


The  amount  of  algebra  involved  in  arriving  at  (56)  is  such  as  to  make  one  susj.ic- 
ious  of  the  result  in  the  absence  of  an  independent  check.  Fortunately,  one  is 
available.  The  case  of  Couette  flow,  solved  by  Serrin  (1959)  and  Joseph  (196o) 
can  be  handled  by  the  same  technique.  The  inner  solution  (29)  remains  the  same, 
as  does  the  outer  solution;  we  have  only  to  set  f  a  1;  the  fact  that  f  varied  af> 
fected  the  matching  only  through  the  value  of  Y.  q.  To  determine  the  value  of 
Y.  o>  have  only  to  determine  the  value  of  y^,  and  here  we  must  replace  our 
previous  condition  (that  i|t  decay  exponentially)  by  a  requirement  that  *  be  sym¬ 
metric  about  the  mid-point  of  the  channel.  Thus  y,.  must  be  fixed  so  that 
3Tt/4-kY_|  or 


3n/4  -  kd/dg”*  -  1)^^^  (y^  -  1)  =  n/2  (57) 

Substituting  in  (54),  we  have 

TT/2-k(l/4B'*.l)‘''^=-(l/4e'*-l)*^^ek/@-(l/4B^-l)*'*^(l/48B^+l)€^k^/2g^+--- 

(58) 

4  1/2 

We  may  thus  go  directly  to  (56),  substituting  n/2  for  3n/4,  and  (1/46  -D 
for  Y.q;  thus,  we  have 

Un\)^  -  (2/tt)(1/4b‘*1)‘^^(1/48B^+1) 

(1/4B^-1)‘''^/T  *  (1/4B^-1)’'^ 

+  (1/4B^-1)‘^^  (1/48B^+1)  ‘^^(3/2^^^)  (n/4)^^^ 

The  same  computations  have  been  carried  out  using  these  equations;  the  results 
are  also  shown  in  Figure  6.  It  may  be  seen  that  the  curve  of  X  vs.  R  extra¬ 
polates  through  the  exact  value  of  R  =  4.  54  as  X-»  0.  while  the  curve  of  k^ 
extrapolates  to  the  exact  value  of  1.  56  at  the  same  value  of  R. 

We  may  conclude  from  this  that  our  estimate  for  X  is  satisfactory  for  X  a  20, 
while  the  k  estimate  is  somewhat  less  reliable,  being  good  for  perhaps  k  a  5. 
The  bending  away  of  both  curves  below  these  points  must  be  attributed  to  pro¬ 
gressive  failure  of  the  approximation. 

DISCUSSION  OF  THE  RESULTS 

Examining  Figure  6,  we  evidently  have  at  R  =  11.  6  a  growth  rate  of  X  ~  51.  5 
and  a  k^  ^  3.  7  (the  latter  value  being  somewhat  less  reliable).  The  critical 
Reynolds  number  for  this  type  of  flow  (corresponding  to  X  =  0)  is  evidently  about 
9.  2,  That  this  should  be  larger  than  that  for  the  Couette  flow  can  be  justiSed  by 
noting  that  in  the  Couette  flow  the  top  halves  of  the  eddies  are  subjected  to  the 
same  strain  rate,  and  hence  are  self-sustaining;  in  the  wall-  layer  flow,  the 
shear  drops  rapidly,  and  the  outer  parts  of  the  eddies  must  be  sustained  by  the 
inner  parts;  hence  the  shear  must  be  proportionately  higher,  or  the  dissipation 
lower,  resulting  in  a  higher  critical  Reynolds  number. 

The  growth  rates  are  quite  large  measured  in  sublayer  variables;  since  sublayer 
time  scales  are  themselves  the  shortest  in  a  wall  flow,  the  evolution  of  these 
eddies  must  be  regarded  as  extremely  rapid.  In  this  connection,  the  observa¬ 
tions  of  IQine  and  his  group  (  e.  g.  -  Kline,  et.  al.  ,  1967)  are  relevant:  wall 
eddies  of  this  general  form  were  observed  to  "burst. "  i.  e.  -  evolve  very  rapidly 
relative  to  other  sublayer  phenomena. 


In  other  reepects  the  dynamics  o£  the  eddies  (from  the  asymptotic  analysis) 
agree  well  with  observationa,  in  particular  the  "mixing  length  approximation" 

(eg.  31a),  which  has  often  been  used  (Townsend.  1956;  Bakewell  and  Lumley, 
1967;  Payne  and  Lumley,  1967). 

The  single  disappointing  feature  of  these  eddies  is  their  size.  Referred  to  sub¬ 
layer  variables.  k4.  =  3,  7  corresponds  to  a  transverse  wave  length  of  roughly 
20,  while  observed  wave  lengths  (e.  g.  Bakewell  and  LunUey,  1967)  are  closer  to 
100.  Our  eddies  have  their  apparent  centers  (the  maxima  of  Ir)  at  about  10, 
while  observed  eddies  seem  to  have  their  centers  closer  to  30.  Evidently  these 
eddies  are  too  small  by  a  factor  between  3  and  5;  can  thi  s  be  explained? 

In  our  considerations,  we  have  not  taken  into  account  the  effect  of  the  presence 
of  the  eddy  on  the  undisturbed  velocity  profile;  to  an  inviscid  first  order  approxi¬ 
mation  the  mean  shear  will  be  reduced  (outside  the  viscous  region  at  the  wall) 
by  an  amount  k‘  /26*  (l/2p^  ~  2  when  R  =  11.  6).  This  reduction  in  shear 
will  favor  the  growth  of  an  eddy  with  the  same  value  of  k^,  but  a  smaller  k  due 
to  the  greater  thickness  of  the  layer.  That  is,  as  the  shear  is  reduced,  all 
length  scales  will  grow  to  keep  R,  kf  roughly  constant.  Hence,  the  Ic  we  have 
found  is  evidently  a  starting  value;  evolution  of  the  eddy  will  involve  progressive 
decrease  in  the  dimensional  k,  until  the  disturbance  energy  reaches  its  peak  and 
begins  to  decay.  Viewed  another  way,  as  the  shear  is  reduced,  the  effective 
Reynolds  number  is  being  reduced*  -  the  eddy  still  feels  the  same  effect  of  vis¬ 
cosity  at  the  wall,  but  has  a  smaller  energy  source;  we  may  expect  the  eddy  to 
move  down  the  curve  toward  X  =  0,  having  at  that  point  a  k+  ~  0.  7  (from  rather 
uncertain  extrapolation  of  the  curve  of  Figure  6).  This  corresponds  roughly  to 
the  value  observed,  in  sublayer  variables,  of  ^  100. 
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APPENDIX  I 


Excerpt  from:  G.  I.  Petrov,  On  the  Growth  of  Qgcillations  in  a  Viscous  Liquid 
and  Treniition  to  Turbulence.  Publication  Number  345,  Trans¬ 
actions  of  the  Central  Aero- Hydrodynamic  Institute,  Moscow, 
1938:  pp.  20-24. 

Lorentz.  selecting  an  initial  velocity  distribution  in  the  form  of  elliptical  vor¬ 
tices,  obtained: 


R  = 


U 

max 


h/v  =  72 


From  a  more  precise  statement  of  the  problem.  Karman  and  Orr  obtained  an 
even  lower  value  Re  ~  42,  which  obviously  does  not  correspond  to  reality. 


Below,  we  will  attempt  to  give  an  explanation  of  this  fact,  and  to  show  how  this 
problem  may  be  correctly  formulated. 


The  linearized  equation  for  the  disturbance  stream  function,  superimposed  on  a 
stationary  plane  parallel  flow,  has  the  form: 


-  R[u  a4*/ax  -  u"aii/dx+  3A*/aT]  =  o.  (8) 

As  we  have  already  shown,  it  is  convenient  to  assume  a  solution  of  the  form 
f(x,y)e*^6*,i,  e.  ,  in  the  form  of  an  oscillation,  the  amplitude  of  which  is  a  func¬ 
tion  of  X  and  y.  The  distribution  of  the  initial  amplitude,  or  the  function  f, 
may  be  obtained  from  the  equation: 

-  R(  u  aAf/ax  -  u"  3f/ax  -  leaf]  =  o  (9) 

and  corresponding  boundary  conditions.  This  equation  is  of  the  fourth  order, 
and  cannot  be  placed  in  self-adjoint  form.  The  coefficients  of  the  equation  con¬ 
tain  two  characterizing  parameters  R  and  0,  of  which  the  first,  in  accordance 
with  its  physical  interpretation,  is  always  real,  while  the  second  is  in  general 
complex.  Thus,  the  problem  may  be  stated  as:  for  a  given  value  of  R  (which 
defines  the  basic  flow)  find  a  system  of  eigenvalues  0  and  a  corresponding  sys¬ 
tem  of  functions  f(x,  y),  assuming  that  an  arbitrary  initial  disturbance  may  be 
expanded  in  a  series  of  these  functions,  or,  in  a  cUfferent  formulation,  find  a 
value  of  R  which  is  associated  in  some  way  with  a  definite  region  of  the  values 
of  0  (for  example,  the  minimum  value  of  R,  below  which  all  given  oscillations 
will  decay,  i.  e. ,  the  imaginary  part  of  0  will  be  negative).  This  value  of  R  is 
ordinarily  used  in  calculating  critical  numbers,  both  as  we  have  seen,  in  con¬ 
nection  with  the  problem  of  stability  of  the  flow  of  a  viscous  fluid  and,  as  it  is 
ordinarily  formulated,  the  determination  of  the  point  of  transition  of  a  laminar 
flow  to  turbulence. 

Neither  of  these  formulations  is,  at  first  glance,  obvious.  The  first,  because 
it  is  not  known  with  certainty  whether  the  complete  solution,  expanded  in  a  ser¬ 
ies  of  periodic  oscillations  of  the  form  fix,  y)  e~'^^  will  decay.  And  the  second 
because  the  critical  number  R  determines  an  initial  growth  diverging  from  the 
real  velocity  of  the  fluid,  and  this  change  in  the  regime  determines  a  change  in 
the  mean  velocity  field. 

Let  us  look  at  several  general  properties  of  our  system.  As  we  have  already  re¬ 
marked,  equation  (9)  cannot  be  placed  in  a  self-adjoint  form,  and  the  eigenv^ues 
of  our  system  will  in  general  be  complex,  in  addition  to  the  obvious  difference 
from  systems  ordinarily  considered  in  the  theory  of  oscillations,  that  the  char¬ 
acteristic  parameters  appear  in  the  coefficients  of  the  derivatives. 

Let  us  examine  the  following  expression: 
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rjT  (x.y)L(£(x.yH<ixdy 


the  integral  being  taken  over  the  region  of  the  flow. 

Here: 

L(f)  =  4if  -  R(U34f/ax  -  u"af/ax  -  ieaf), 

i.  e.  ,  the  right-hand  («ic)  aide  of  equation  (9). 

Thia  expreaaion  can  be  tranaformed  to  the  following  form; 

J‘JT(x.y)L(f{x,y)dxdy  a  JJafAfdxdy  +  Rj*]"  (U  af/ax  Af  + 
U  T  hf/Zx  +  igf  Af)<lxdy  +  •  •  • 


(20) 


plua  a  contour  integral. 

If  the  function  f(x.y)  satiafiea  equation  (9),  then  expreaaion  (20)  muat  equal 
zero,  and  we  may  determine  R  aa  the  ratio  of  the  two  integral  expreaaions, 
and  aeek  the  critical  value  aa  the  minimum  of  thia  ratio,  aetting  gj  equal  to 
aeveral  valuea. 

Since  we  have  given  boundary  confitiona,  let  ua  examine  the  variation  of  juat  the 
double  integrala.  The  problem  of  the  determination  of  the  critical  value  may  be 
aet  aa  an  iaoperimetric  problem,  i.  e.  ,  to  find  the  extremum  of  the  function^ 

If  s  JJhrAfdxdy  under  the  condition 

I2  *  /J  •  [  U  sT/ax  Af  +  U  f  af/ax  +  Igf  Afjdxdy  =  1. 

The  Euler  equation  for  our  problem  will  be  equation  (9),  and  for  f  the  equation 
of  ita  adjoint; 

AdT  -  X  [u  aA7/ax+  2u'  a^7/axay  -  igAlj  =  o.  (21) 


But,  in  order  that  our  integrala  have  the  aame  phyaical  interpretation  aa  in  the 
work  of  Lorentz,  Karman,  et.  al.  which  we  deacribed  above,  we  muat  aasume 
that  T  ia  the  complex  conjugate  of  a  complex  function  f,  i.  e.  ,  that 


f  =  f  +  if. ,  to  7  =  f  -  if. . 
r  1  r  i 

All  of  our  concluaiona  follow  from  the  form  of  the  ititial  diaturbance,  i.  e.  ,  by 
taking  into  account  that  the  coefficjenta  of  our  equation  do  not  depend  on  x. 
taking  f  =  0e'®*  and  then  f=J(y)e'*®*(where  7(7)  i*  fhe  complex  conjugate  of  0(y)), 


Then  the  expreaaion  (2)  takea  the  following  form: 

qJ'  (5’"0"+2a^0'0'+a*00) <ly+Rj’[aU  l/2(00'-00')  -g.(@'0'+a00)  J  dy 

-  iaR  jjJ*  [(U  -  c^)  (0'0'  +  a^00)  +  1/2  U"00]  dy 


(20a) 


The  integrala  are  taken  over  a  region  contained  bwtween  the  walla,  and  over  a 
wavelength  of  the  given  diaturbance,  ao  that  the  contour  integral  will  vaniah. 


The  real  part  of  expreasion  (20a)  can  be  regarded  as  the  change  in  energy  in  a 
strip  of  one  wavelength  of  the  disturbance,  and  may  be  obtained  by  examining 
the  expression  for  dE/dt  dx  dy  (see  expression  (6)).  taking  into  account  that  the 
real  part  of  the  stream  function  is  4  =  0  (y)e^((ix-gt)  ^nd  integrating  in  x  over 
one  wave  length. 

The  same  functional  can  be  obtained  as  an  expression  (^L(0)  +  0L(0)  )dy, 
where  L(0)  is  the  left-hand  side  of  equation  (13),  while  C  is  the  corres¬ 
ponding  expression  for  the  conjugate  function,  and  has  the  form: 


L  (0)*  0 


2o^0"  +  a^0  +  iaRt(tJ  -  c! 


)  (0"  -  -  U"0], 


i.  e.  -  the  left-hand  side  of  equation  (13)  with  the  complex  parameters  changed 
to  their  conjugates. 

The  imaginary  part  of  expression  (2)  may  be  obtained  as  an  expression 

(0  L  (0)  -  0E  (0) )  dy 

Both  of  these  expressions  for  functions,  satisfying  equation  (13),  must  vanish. 
Then  from  the  real  part  we  may  express  R,  or  0  for  given  R,  as  a  ratio  of 
functionals.  The  imaginary  part  satisfies  a  Zol  'berg  condition,  from  which  it 
follows  that  for  undisturbed  velocity  profiles  not  having  an  inflection  point,  c^ 
is  always  less  than  U^iax'  may  obtain  the  phase  velocity  c^  or  the  wave 

number  9;  as  a  ratio  of  functionals. 

But  in  order  that  the  disturbance  be  dynamically  possible,  i.  e.  ,  satisfy  equation 
(13),  the  real  and  imaginary  parts  must  vanish  simultaneously,  and  the  problem 
of  the  critical  yalue  must  be  formulated  as  a  problem  of  the  extremum  of 
oj*  nr'0"  +  2a^7'0'  +  a’T<S)<iy  under  the  condition 

qJ*^{1/2  aU(0'0- 00')  +  iai  (U-c^)(0’0'+a^00)+l/2  U  0'0  |}dy=l 

(which  includes  the  condition  that  the  imaginary  part  vanish).  The  Euler  equa¬ 
tions  for  this  problem  are  equation  (13)  for  0  and  its  adjoint  for  0,i.  e.  , 

—  IV  2—  4—  —  2—  — 

0  -  2a  0"  +  0  0  -  iaRt(U  -  c)  (0"  -  a  0)  +  2U'0’]  =0. 

But,  if  0  satisfies  equation  (13),  then  the  function  must  satisfy  the  equation 
‘E(0)  =  0.  This  is  possible  only  if  the  following  relation  holds: 

(U  -  c^)  (0"  -  a^)  -  1/2  U"0  +  U'0'  =  0,  (22) 

i.  e.  -  if  equation  (13)  and  its  conjugate  are  anti- symmetric. 

Bat  the  condition  (22)  very  much  reduces  the  class  of  functions  which  can  give 
our  solution.  Since  the  solution  of  equation  (22)  depends  only  on  two  arbitrary 
parameters,  the  solution  of  the  variational  problem  in  such  a  formulation  is  ap¬ 
parently  impossible.  This  explains  the  failure  of  all  the  attemps  to  resolve  the 
problem  of  the  so-called  energy  method,  since  apparently  disturbances  obtained 
from  the  minimum  value  of  the  ratio  which  determines  R,  are  dynamically  im¬ 
possible. 

Thus,  the  critical  number  R  may  be  found  as  the  smallest  value  of  the  ratio 
Il/R(l2)  (where  R(l2)  is  the  real  part  of  the  complex  quantity  I2),  in  which  are 
substituted  solutions  of  equation  (9),  but  these  functions  do  not  give  an  extremum 
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of  this  ratio.  It  may  aiao  be  said  regarding  the  eigenvalues  0  =  0^  +  i0.,  that  in 
the  system  under  consideration  they  apparently  do  not  possess  an  extremal  char' 
acter. 


Figure  1  Reichardt's  (1951)  Viscosity  Expression. 


Figure  3  Structure  of  ^  and  Corresponding  Cells  for  One  and  Two  Turning  Points.  Note  Change  due  to  Viscosity 
(dashed  curves). 
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